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Abstract: The persistent discrepancy between observations of "^Li with putative primor- 
dial origin and its abundance prediction in Big Bang Nucleosynthesis (BBN) has become a 
challenge for the standard cosmological and astrophysical picture. We point out that the 
decay of GeV-scale metastable particles X may significantly reduce the BBN value down to 
a level at which it is reconciled with observations. The most efficient reduction occurs when 
the decay happens to charged pions and kaons, followed by their charge exchange reactions 
with protons. Similarly, if X decays to muons, secondary electron antineutrinos produce 
a similar effect. We consider the viability of these mechanisms in different classes of new 
GeV-scale sectors, and find that several minimal extensions of the Standard Model with 
metastable vector and/or scalar particles are capable of solving the cosmological lithium 
problem. Such light states can be a key to the explanation of recent cosmic ray anomalies 
and can be searched for in a variety of high-intensity medium-energy experiments. 
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1. Introduction 

Rapid progress in observational cosmology during the last decade brought about the mea- 
surements of many cosmological parameters, including a precise determination of the 
baryon-to-photon ratio r]b [I] by the WMAP satellite experiment. This puts the pre- 
dictions of the Big Bang Nucleosynthesis (BBN) theory on firm footing and allows for less 
ambiguous comparison with observations. The current status of standard BBN (SBBN) 
with the input from WMAP can be summarized as follows: there is no dramatic (order 
of magnitudes) disagreement between predictions and observations, but there is no ideal 
concordance either. For a series of recent reviews see e.g. [2, 3, 4, 5]. 

There are at least two quantitative problems that look worrisome: different mea- 
surements of the deuterium abundance, although on average consistent with the SBBN 
prediction, exhibit a significant scatter. This scatter may be the sign of underestimated 
systematic errors or the manifestation of significant astration, thus hinting on a higher 
primordial value for the deuterium abundance. In contrast to deuterium, the scatter of 
^Li/H data points along the so-called Spite plateau [6] is rather small, which for a long 
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time thought to be a compelUng argument for the primordial origin of Li in these obser- 
vations. As is weh-known, this value is a factor of 3-5 lower than the SBBN prediction, 
'''Li/H = (5.24lo;g7) x 10"^° [3], which is the essence of the lithium problem. 

How serious are these problems of SBBN? It is entirely possible at this point that 
future highest quality observations of D/H in quasar absorption clouds would render D/H 
in accordance with SBBN in combination with less scatter. Moreover, more elaborate stellar 
evolution models with ab-initio calculations of lithium diffusion may point to a systematic 
and uniform reduction of the SBBN value by a factor of 3 or so. At this point, it is too 
early to declare SBBN being in serious trouble. However, it is also tempting to speculate 
that some subtle particle physics interference in the early Universe may have resulted 
in distorted abundances of the primordial elements, and possibly led to the reduction of 
lithium abundance. 

The primordial value of the lithium abundance is given by the freeze-out BBN value of 
^Be-|-^Li, with atomic ''Be decaying to lithium at the later stages of cosmological evolution. 
The current lithium problem stems from the overproduction of ^Be at T ~ 40 keV. '^Be 
cannot be destroyed by protons directly, but instead is depleted via the following chain of 
exoergic reactions 

^Be + n ^ '^Li ^ ''He + %e. (1.1) 

At the second step ''Li is destroyed by proton reactions, which remain faster than the 
Hubble rate until T ~ 10 keV. 

Different classes of modified BBN models where an additional reduction of lithium can 
happen were analyzed in the literature over the years. As was first pointed out by Reno 
and Scckcl [7], a (moderate) injection of "extra neutrons" around the time of formation of 
''Li and ^Be leads to an overall depletion of ^Li-l-'^Be by intensifying the destructive chain 
(1.1). This was emphasized again after the CMB determination of r]b in Ref. [8], where it 
was demonstrated that any particle physics source is capable of reducing ^Li-|-^Be as long 
it leads to the injection of O(10~^) neutrons per nucleon at relevant temperatures. Perhaps 
the most natural source for a neutron excess is the decay of massive particle species X. 
Independently from the motivation of reducing the ^Li abundance, a lot of work has been 
invested into BBN models with decaying or annihilating particles releasing a significant 
amount of energy into the primordial plasma [9, 10, 11, 12, 13, 14, 15]. 

So far, most of the analyses have concentrated on the injection of energy by relics 
with masses comparable to the electroweak scale. This is largely motivated by theoretical 
arguments in favor of new physics at and below the TeV-scale, and by the possibility of 
having dark matter in the form of weakly interacting massive particles (WIMPs). It can 
be easily shown that it is unlikely that the residual WIMP annihilation is responsible for 
the reduction of ^Be-|-^Li, simply because the total energy injected via such mechanism is 
way below the required levels [7, 16]. Known examples that "work", i.e. scenarios in which 
the ^Li abundance is reduced while other elements are still agreeing with observations, 
typically deal with unstable weak-scale particles. These include some supersymmetric 
scenarios with the delayed decays of charged sleptons to gravitinos [17, 13, 18]. The source 
of extra neutrons in these models is linked to the presence of nucleons among the decay 
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products. An alternative plausible mechanism for reducing ^Li is the catalytic suppression 
of the ^Be abundance by the capture of massive negatively charged particles [19, 20, 21, 22], 
that are again linked to the weak scale. In this paper we investigate whether the suppression 
of the lithium abundance can be triggered by the decays of metastable GeV and sub-GeV 
scale electrically-neutral particles. We address two types of models: the WIMP-type where 
particle decays were preceded by the depletion through the annihilation, and the super- 
WIMP type, where the abundance of decaying particles is set by the thermal leakage of 
Standard Model (SM) states into an initially vacuous super- WIMP sector. 

During the last two years, the GeV-sector phenomenology experienced some degree of 
revival due to its possible connection to the enhancement of the leptonic fraction of cosmic 
rays in the pair annihilation of dark matter WIMPs [23, 24]. Particularly noteworthy arc 
the enhancement of the positron fraction seen by the PAMELA satellite experiment [25] 
and the harder-than-expected spectrum for electrons and positrons detected by the FGST 
instrument [26]. In such scenarios the GeV-scale particles are designated as "mediators" 
connecting the dark matter and Standard Model sectors [27], allowing to seclude the dark 
matter by choosing the the SM-mediator coupling to be very small, but at the same time 
keeping the galactic annihilation rate enhanced over the "standard" WIMP scenario. 

The motivation for GeV and sub-GeV scale mediators comes from the following con- 
sideration: the lightness of a gauge boson mediating an attractive force in the dark sector, 
V{r) = —an/r x exp(— myr), compared to the characteristic WIMP momentum inside 
the galaxy, moVgaX, ensures a Coulomb/ Sommerfeld enhancement of the annihilation cross 
section relative to its freeze-out value. Choosing my <mDV^ 100 MeV ^ 1 GeV, results 
in a 'Kao/v enhancement of the cross sections at the relevant velocities. Moreover, once 
recombination to WIMP-bound states becomes kinematically possible, my < rriDaj^/A, 
this process dominates numerically over the direct two-body annihilation. The resulting 
annihilation cross section for fermionic WIMPs can be enhanced over the freeze-out value 
by two-to-three orders of magnitude [24] : 



{'^v)f.o. Vgal 



~ 0(10^ - 10^). (1.2) 



The numerical enhancement due to the bound-state effect over the naive Sommerfeld value 
is about a factor of 7, accompanied by a possible additional enhancement due to an in- 
creased lepton multiplicity in the final state with angular momentum J = 1. We note 
in passing that the important effect of WIMP-onium formation was missed in a recent 
re-analysis of Ref. [28] , that led to the erroneous conclusion that the enhancement factor 
remains smaller than 100 over the whole parameter space; larger factors were found in [29]. 

Is it reasonable to expect that the same models that fit the PAMELA and FGST signals 
[23, 24] are also responsible for the suppression of ^Li abundance? Even with the inclusion 
of some very generous enhancement factors, the energy injection during BBN triggered 
by the annihilation of electroweak-scale WIMPs remains rather small. Moreover, models 
designed to explain the PAMELA signal tend to minimize the fraction of baryons/anti- 
baryons in the final state [23, 24, 30, 31]. All that, together with previous investigations of 
BBN with annihilation-induced energy injection, tend to indicate that WIMP annihilation 
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itself cannot be used as a solution to the lithium problem. Therefore, the only chance of 
altering the BBN predictions for ^Be+^Li within this class of models is if the GeV-scale 
particles from the mediator sector are themselves relatively long-lived, and their decays 
happen during or after BBN. Interestingly enough, it turns out that the minimal ways of 
coupling the Standard Model to light mediators often implies the longevity of particles in 
the GeV sector [32]. 

The main mechanism by which the decays of the GeV-scale relics in the early Universe 
can influence the outcome of the BBN is the injection of light mesons such as pions and 
kaons as well as muons that all lead to the extra source of p — > n conversion. In this 
work we explore such scenarios, finding the "required" number of injected Tr~ , K~ , and 
fjr triggering p ^ n conversion in the right amount, as well as the "optimal" lifetime- 
abundance window for such injection. In order to have a consistent cosmological picture, 
the abundance of parent GeV-scale relics should be small enough so that they carry only a 
small fraction of the energy density of Universe during BBN, but still provide a noticeable 
number of mesons and muons per nucleon. We show that, in fact, many models with 
GeV-scale relics fulfill this requirement, including some variants of the models designed 
to fit PAMELA and FGST signals. We also find that both the WIMP and super- WIMP 
modifications of a secluded U(l)5-model is capable of reducing the lithium abundance to 
the observable level. 

The structure of this paper is as follows. The next section contains the analysis of the 
injection of tt, K, and particles vs. the timing of injection that is needed for reducing ^Li. 
Section 3 investigates a variety of different models for decaying GeV and sub-GeV scale 
relics, and finds the regions of parameter space that lead to the depletion of the lithium 
abundance. Our conclusions are summarized in Sec. 4. Appendix A contains relevant 
details with physics input that went into our BBN code. 

2. Meson and neutrino injection during BBN 

We begin this section by presenting an overview of the physics processes triggered by the 
decays of meta-stable GeV-scale relics X during BBN. Thereby we shall estimate timescales 
and interaction efficiencies of the crucial reactions which eventually lead to the reduction 
of the overall ^Li BBN prediction. Subsequently, the various final states in the decay of X 
are considered in detail. 

The central parameter entering the discussion is the Hubble expansion rate H as it 
normalizes any interaction rate during nucleosynthesis. The most relevant epoch for BBN 
corresponds to the time bracket of 100 sec < t < 1000 sec or, equivalently, to a temperature 
range of 



where Tq is the photon temperature in units of 10^ K. Thus, in the interesting regime after 
the annihilation of the electron-positron pairs, the Hubble rate is given by 



1.2 < Tg < 0.4, 



(2.1) 





(2.2) 
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where gcs — 3.36 counts the radiation degrees of freedom; Mp ~ 2.43 x 10^^ GeV is the 
reduced Planck scale. 

The synthesis of ^Be, proceeding via ^He + ^He — > ^Be + 7, occurs in a more narrow 
temperature interval, Tg c± 0.8 0.4, in which the Hubble rate is ~ 10~^ sec~^. This 
occurs immediately after the opening of the deuterium "bottleneck" at Tg ~ 0.8, when 
significant quantities of ^He and ^He are formed, with a helium mass fraction of ~ 0.25 
and a number density of ^He per proton of YsHe — 10~^. The rate for ^Be formation per 
alpha particle, 

r.Be - (10-=^ sec-^) Tl/'e-^'-'/^^'' {Tg < 0.8) (2.3) 

always remains slower than the Hubble rate and quickly becomes completely inefficient 
due to a strong exponential Coulomb suppression. On the other hand, a non-standard 
"thermal" neutron abundance, at the level comparable to lO^'^ has the opportunity to 
reprocess some fraction of ^Be via ^Be + n — > ^Li + p since the rate for a neutron capture 
per ^Be nucleus is 

rrBe->rLi ^ (0.4 sec"^) T| . (2.4) 

Indeed, looking at the fraction /rgg of ''Be which can, in principle, be converted within a 
Hubble time Atn, 

/.Be -J-^AtH -H-'r,^,_ru ~ -W'TgY^, (2.5) 

shows that already ~ 5 x 10~^ at Tg ~ 0.5 can induce an 0(l)-conversion of ''Be. ^Li 
in the final state of this reaction is then quickly burned by protons. 

The neutron participation in the lithium depleting chain does not influence the abun- 
dance of neutrons themselves, because of the small abundance of ^Be. Instead, "extra 
neutrons" are depleted by protons via the p + n — )■ D + 7 reaction, which, at the rele- 
vant temperatures, remains faster than the rate for neutron decay. Comparing the relative 
changes of ''Be and deuterium, 

fr^ = _^ r7B,^^Li ^ _^q5 p/jj) (2.6) 

one can infer that 0(l)-depletion of ^Be due to an excess of neutrons will be accompanied 
by an 0(l)-rise in the deuterium abundance since D/H ~ 10^^. Therefore, the viability of 
the "extra neutron" solution to the lithium problem should be judged form more accurate 
quantitative investigations of such scenarios. 

The decay of GeV and sub-GeV relic particles with a rate comparable to (2.2) at 
Tg ~ 0.5 will not lead to a significant population of energetic photons and electrons. The 
reason for that is well-known: very abundant and energetic photons and electrons degrade 
the energy of the decay products well below the nuclear dissociation thresholds (see e.g. 
[10]). Since we explicitly assume the absence of nucleons and anti-nucleons in the final state. 
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the main effect on the freeze-out abundances of hght elements will come from injection of 
mesons and neutrinos, that are capable of triggering the p ^ n conversion and reducing 
^Li+^Be. In the following subsections we consider these mechanisms in turn, starting from 
the simple estimates of the required number of mesons and neutrinos. 

2.1 Estimates on meson injection 

Among the the decay products of the GeV-scale metastable states only the long-lived 
mesons and neutrinos will have a chance to interact with the light elements. Among those, 
the most important are charged pions, tt^, as well as charged and long-lived kaons, 
and Kl, respectively. Their masses and lifetimes are given by [33] 





= 139.6 MeV, 


r^± 


= 2.60 X 10"^ 


sec. 


(2.7) 


mK± 


= 493.7 MeV, 




= 1.24 X 10"^ 


sec. 


(2.8) 




= 497.6 MeV, 




= 5.12 X 10"^ 


sec. 


(2.9) 



The decay rate of meson i is related to its lifetime at rest by a familiar time-dilatation 
formula: 

rLc = (^)^^fewxl0^sec-i (2.10) 

where in the last equality we have taken 7 ~ 1. More accurately, one has to consider the 
average value {rrii/Ei) over the "lifetime" trajectory of injected mesons. The averaging 
procedure is significantly different from the usual thermal average when particles obey 
a Maxwell-Boltzmann distribution. It is determined by the efficiency of various energy 
degradation mechanisms, of which the most important are Coulomb scattering on electron 
and positrons and inverse-Compton-type scattering on background photons. Together they 
determine the thcrmalization (or stopping) rate Fg^^p. The stopping rate may have a direct 
impact on the the reaction rates with nuclei F^y^, which are of the primary interest in this 
paper. 

If injected sufficiently early, the charged mesons are effectively stopped within time 
intervals shorter than their lifetimes Tj, and before having a chance of participating in a 
nuclear reaction. The comparison of stopping and decay rate is given by the ratio 

"^""^ - '^^^ - frj,,, r -1^ V\ (2.11) 



where rgtop measures the energy degradation time from injection energy £^0 to some energy 
-Emin, below which the kinetic energy of mesons is irrelevant. dEi/dt represents the energy 
degradation rate of the particle i traversing the BEN plasma. Requiring complete stopping 
would correspond to the choice of £?min ~ 3r/2, but a more relevant parameter is some 
characteristic nuclear energy scale. For reactions leading to dissociation of ^He, .Bmin 
would correspond to -Emin ~ 20 McV. In the temperature window (2.1) and for injection 
energies of E'o ^ 1 GeV, Coulomb interactions on give the dominant contribution to 
dEi/dt (for details see e.g. [10]). Since the number densities of electrons and positrons 
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drop exponentially, the inverse Compton scattering on background photons becomes the 
dominant energy loss mechanism for T < 25 keV. Taking into account both channels for 
energy loss, one typically finds that 

^ < 1 for Tg > (0.3 H- 0.4). (2.12) 



T 

stop 

This point is illustrated in Fig. 1, where the dividing line of "1" separates the two regimes of 
complete and incomplete stopping. The exponential sensitivity to temperature is reflected 
in almost vertical contours. Figure 1 tells us that for t < 10^ sec the charged mesons are 
thermalized before they decay. Thus, for i = 7r='= (and analogously for K"^) injected around 
Tg ~ 0.5 the following hierarchy of scales is applicable: 

ii«r-«r3,,<r-„p. (2.13) 

The rate for charge exchange reactions with the most abundant nuclear species, i.e. 
the rate of proton-to-neutron conversion, is given by 

r; = n,{av);^ (3 x 10^ sec'^) (2.14) 

The averaging procedure is again determined along the trajectory of injected mesons and 
includes both the "in-flight" and "at-rest" contributions. However, for the purpose of a sim- 
ple estimate, Eq. (2.14) uses the normalization on a typical size of a pion induced reaction 
at the threshold. In case of incomplete stopping, this will underestimate the proton-neutron 
conversion because of the delta-resonance enhancement of the charge-exchange reaction. 
Equation (2.14) immediately leads us to the probability of charge exchange reaction of a 
stopped pion (kaon) on protons during the meson lifetime, per each injected tt~{K~) at 
Tg ~ 0.5: 

/•oo 

Pp%„ = / exp(-rdec(t - iinj))rpdt ~ t;t,± ~ 0(10-^) (2.15) 

Since every charge exchange reaction leads to the creation of an extra neutron, the 
estimate (2.15) represents the efficiency of producing neutrons from negatively charged 
kaons and pions. Notwithstanding a rather crude nature of this estimate, the probability 
(2.15) tells us that an injection of 0(10) negatively charged pions per proton should be 
equivalent to the injection of 0(10"^) extra neutrons, and thus capable of reducing the 
overall ^Li-|-''Be abundance by an 0(1) factor. 

Since for kaons the charge exchange cross sections are even higher [7], one can achieve 
an adequate suppression of ^Li-|-^Be even with the injection of one K~ per baryon. The 
case of X-decays into neutral kaons is special with respect to other meson final states 
in the sense that has a relatively long lifetime (2.9) but, unlike tt^ or is not 
stopped by electromagnetic interactions. From the conservation of isospin and from charge 
independence of strong interactions, we nevertheless expect the impact of Kl on the BBN 
predictions to be similar to that of K~ . 
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Figure 1: Stopping of cliarged pions and kaons in the plasma as a function of temperature Tg, 
initial injection energy Ei^, and the minimal energy below which the kinetic energy of mesons can 
be neglected, E'low The contour Tdoc/Tstop = 1 is the dividing line between efficient (left) and 
incomplete (right) stopping. 



In our treatment, the extra meson species are included in the set of Boltzmann equa- 
tions, along with the population of parent particles X that decay into mesons and have an 
abundance with the simple exponential time-dependence, 

Yx{t) = Y^eM-t/Tx), (2.16) 

with Tx being the X-lifetime. We choose to normalize the X-abundance on the total 
number of baryons, Yx = nx/ni^; in the following we shall often drop the superscript 
on for simplicity. When the fast stopping of charged particles is operative at early 
times, the amount of the injected energy and therefore the actual mass of the X-relic 
are not entering the problem (with X — t- K^K^ being the exception.) For later times, the 
incomplete stopping brings the dependence on the injected energy. We write the Boltzmann 
equation for meson species i in the following form 

e?V,rf,j - y.rji,, - y^r^. (2.17) 

j 

where Tf^- is the rate of i-injection from source j with multiplicity When considering 
only the primary meson production from X-decays such as in the case of , Ff^- ~ 
= 1/tx- Given the hierarchy of interactions (2.13), the last term in (2.17) is small. 
Both r^g^ and include the effects of incomplete stopping that depends on the initial 
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energy injection Eq and background temperature at the time of injection, 



ri^, = r^,SEo,T); r% = T%{Eo,T). (2.18) 

Neglecting subtleties of incomplete stopping at this point, one can find an approximate solu- 
tion to the Boltzmann equation for mesons in the quasi-static approximation dY^-^^/dt = 0, 

yr = E d^V,r/„./r^,, ~ 0(10-10) X Yx, (2.19) 
j 

where in the second relation we assumed multiplicities to be order one, and took Finj ~ 
10"^ sec. Even though the equilibrium values for the meson abundances remain low at all 
times, the p — >■ n conversion rate for Yx ~ 0(10) is only five orders of magnitude slower 
than the Hubble rate, thus causing 0(10"^) protons be converted to neutrons within one 
Hubble time. This is consistent with our estimate of efficiency (2.15). 

Finally, the pion(kaon)-induced transitions ^Be — > ^Li and ^He T also lead to 
the depletion of ^Be. However, these processes are far less important compared to p — )■ 
n interconversion. If we assume a fractional O(10^^)-conversion of p — )■ n due to vr" 
charge exchange, a similar figure would stand for the ''Be — )■ ^Li and ^He — )• T conversion 
probabilities, which is completely negligible in the final lithium abundance. 

2.2 Estimates on muon and neutrino injection 

In order to account for the injection of muons and neutrinos one has to deal with quite 
different physics. The charge exchange reactions of muons on nucleons arc mediated by 
weak interactions and can be neglected during the muon lifetime of two microseconds. 
However, the decays of muons, fi~ — > u^v'(.e~ , source energetic electron antineutrinos i/g, 
that survive for a long time, which increases their probability of charge-exchange interaction 
with protons. 

To make our discussion more concrete, we shall assume that all neutrinos originate 
from muons decaying at rest so that Ei, < m^/2. In this case, the relation between the 
relevant rates is quite different from (2.13): 

r;, r,-top « H, (2.20) 

where Tg^^p is the rate of antineutrino energy degradation due to scattering on background 
neutrinos and electron-positron pairs. This rate scales with temperature and the energy of 
non-thermal neutrinos E^iT) injected at Tinj as follows: 

r^top - GlE^{T)T' ^ GlT'^-^^p^, (2.21) 

-tmj 

where we disregard the difference between neutrino and photon temperatures. It is indeed 
much smaller than the Hubble expansion rate. 
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where 3 MeV scale enters as the decouphng temperature of background neutrinos, and the 
energy of non-thermal antineutrinos is taken to be the average energy in the muon decay, 
(^.(Tinj)) = ^m^ ~ 32 MeV. 

The charged current rate for antineutrino interactions with protons is given by 

rl = ^ lO-"' cm^ X ^jj^ll^ (3.6 X 10-" sec-') x (2.23) 

where Ei, mg is assumed. The ratio of to the Hubble rate gives the efficiency of 
producing extra neutrons from each neutrino injected with 0(30 MeV) energy: 

This efficiency is several orders of magnitude smaller than Eq. (2.15), and we conclude that 
O(lO^) muon decays per proton are required around Tg ~ 0.5 in order to produce O(l0~^) 
extra neutrons. At the same time, Pp^n ^ -^pUn introduces a natural simplification to the 

problem, as the effects of secondary neutrinos originating from pion and kaon decays can 
be ignored relative to the direct influence of vr^ and via the charge exchange reactions. 
Moreover, as the efficiency (2.24) also suggests, a direct nuclear-chemical impact of electron 
antineutrinos on the lithium abundance via ^Be + i>e ^ ^Li + e'^ {Q = —0.16 MeV) is 
negligible. 

Having checked that the efficiency of neutron production from pions, kaons and muons 
(neutrinos) can be sufficient for the resolution of the lithium problem, we now turn to 
more detailed calculations. Some further details on the Boltzmann code that we use can 
be found in the Appendix. 

2.3 Decays to pions 

Charged pions are a likely final state for almost any hadronic final decay mode. For 
example, even if X would exclusively decay into K^, i.e. = in Eq. (2.17), pions 

would still be populated by subsequent K='=-decays, Yp-± 2± Y^±, since ^ = 0{l). 

Charged pions have a chance to interact with protons and ^He before decaying. Fully 
thermalized pions have the following charge exchange reactions on nucleons with positive 
energy release: 

TT'+p^n + j: {av);~^^^ 0.57 mh, Q = 138.3 MeV, (2.25) 
7r-+p^n + 7r°: (cti>);^„(^o) ^ 0.88 mb, Q = 3.3 MeV, (2.26) 
TT^ + p + TT^ : {av)l~ 1.7 mh, Q = 5.9 MeV. (2.27) 

These reactions interconvert neutrons and protons and thereby increase the n/p ratio, 
because protons are far more numerous once neutrons have been incorporated into ^He. 
The threshold value for the +p reaction cross section can be extracted from the lifetime 
of the pionic hydrogen atom (see e.g. the review [34]). The strength of the channels tt^ and 
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7 at threshold can be inferred from the Panofsky ratio [35]. Our value at the (thermal) 
threshold 

' <^<nM + i<^<n(.")] ^ Fv^- X 1-45 mb, (2.28) 



is in good agreement with the one used in the BBN literature [7] from where we also 
took the cross section for vr^ + n. This threshold value is subject to the usual Coulomb 
enhancement that is accounted for with the multiplicative factor Fp^- in (2.28). For the 
negatively charged pions and at the temperature range of interest (Tg ~ 0.5) this translates 
into an enhancement of the reactions rates by a factor of -F],^- — 2. Further details on 
Coulomb enhancement are provided in the Appendix A.l. 

The incomplete stopping of pions introduces an important modification to the rates 
of these reactions, and to the efficiency of p — )• n conversion. To quantify this effect, we 
introduce the correction factor k{Eq,T) as a function of primary kinetic injection energy 
Eq that includes both the effect of the time dilatation (and effectively longer lifetimes of 
fast pions) and, more importantly, the momentum dependence of the cross section that has 
a strong peak around a pion energy of ^ 180 MeV: 

.(^„, . ^ f exp (- ^f,^) . (-9) 

The exponential factor in this expression is the survival probability, that is the probability 
of finding a meson at time t after its injection at t = 0. Defined this way, the correction 
factor k{Eq,T) = 1 in the limit of Tstop ^ ^tt- Notice also that due to the scale hierarchy, 
H <^ , we can disregard the effects related to cosmological expansion inside the integrals 
of Eq. (2.29), and set the initial moment of injection to be t = 0. 

The kinetic energy E^^ of an injected pion, and the time t along the "lifetime" trajectory 
are related via the rate of energy loss: 

/•^o dE' dE (dE\ fdE\ 

As alluded before, dE/dt is a strong function of the background temperature. Finding the 
explicit dependence between time and kinetic energy of injected pions numerically, and 
using the experimental data for the inelastic p + 7r~ reaction away from threshold, we find 
the correction factor for the efficiency, k{Eq,T). Figures 2 and 3 plots this factor for the 
representative temperature and initial energy slices. One can see that due to incomplete 
stopping at late times, the correction factor can be very large, possibly reaching 



(o-^')th 

where 7max — 2.3 corresponds to the pion momentum at delta-resonance, where the ratio 
of cross sections reaches (o"v)max/ (c"i')th — 14. 

At the next step we account for the possible reactions on ^Hc. Again, we separate our 
discussion into the reactions induced by stopped and in-flight pions. Charge exchange re- 
actions of thermal tt^ on ^He are not possible because of the deep binding of '^He compared 
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Figure 2: The correction factor k{Eq,T), relating the efficiency oi p ^ n conversion by pions 
in-fiight to the thermal function of temperature for representative values of injection 

energy Eq; Eth = 3T/2. The dashed line shows the case of maximal efficiency due to pion injection 
at the delta-resonance energy. Formally, k < 1 is possible for low temperatures and small injection 
energies because the Coulomb corrections of a fully thermalized pion are substantial. 

to other mass-4 isomers. However, vr^ can be fully absorbed by ^He leading to ~ 100 MeV 
energy release and a variety of nuclear final states. Moreover, the reactions with thermal 
TT"*" are very suppressed because of the Coulomb repulsion, and thus we concentrate on vr" 
absorption. A measurement of the ground state level width Tis = (45 it 3) eV [36] of pionic 
helium then allows us to obtain the total low-energy in-flight cross section (crv) ~ 7.3 mb; 
for further details see Appendix A. 2. The branching ratios into the different final states 
have been measured [37] to be (Tn) : (Dnn) : (pnnn) = (17±9)% : (63±26)% : (21±16)%. 
Adopting the central values (with 20% pnnn) we arrive at 

7r-+^He^T + n : (o-v)f.„ ~ 1.1 mb, Q = 118.5 MeV, (2.32) 
7r"+^He^D + 2n: (o-t;)5„„ ~ 4.1 mb, Q = 112.2 MeV, (2.33) 
7r-+*^He^p + 3n : {av)^~^^ :^ 1.3 mh, Q = 110 MeV. (2.34) 

Before using (2.32) in our code, we account for Coulomb corrections, which leads to the 
enhancement of F4jje7r- — 3.5. 

The reactions of "in-flight" pions with ^He can also be significantly enhanced by the 
presence of the delta-resonance. Again, one has to distinguish two types of reactions, 
inelastic scattering, vr^ + ^He — t- vr + A^, and absorption, vr^ + '^He — t- N where N represents 
a variety of multi-nucleon/nuclear final states. We use the results of experimental studies 
[38, 39, 40] to account for the pion-^He reactions across the delta-resonances. We also 
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Figure 3: The correction factor k{Eq,T), relating the efficiency oi p ^ n conversion by pions 
in-flight to the thermal function of injection energy Eq for Representative values of 

temperature. The curves corresponding to late injection Tg = 0.15 and 0.3 when stopping does not 
occur follow the familiar energy dependent profile of the pion-nucleon cross section with the broad 
delta resonance at 180 MeV. 



extrapolate these results to the threshold regions, and calculate effective cross sections 
to the various (exclusive) final states by averaging over the pion "lifetime" trajectory 
while taking into account the respective energy-dependent branchings. The details of this 
procedure are summarized in Appendix A. 3. 

Besides enhancing the number of free neutrons, the reactions on '^He have an additional 
important effect: they produce energetic ^ = 3 elements that are able to participate in the 
non-thermal reactions leading to ^Li. For example, reactions with thermalized vr" in Eq. 
(2.32) in ~ 17% of all cases contain T nuclei with an energy E^^ ~ 30 MeV injected into 
the plasma. This leads to a possible secondary source of ^Li via their fusion on ambient 
alpha particles [7, 41], 

T + ^He|bg ^ ^Li + n Q = -4.8 MeV. (2.35) 

The number of produced ^Li per injected T (and likewise per injected ^He) can be found 
by tracking the T-degradation from i?T,in until the threshold energy in the frame of the 
thermal bath, E]'^. ~ 8.4 MeV, 

where vt is the velocity of T. As is well-known, the production of ^Li becomes more 
efficient at late times, when its thermal destruction slows down. This occurs in the regime 
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of <C Tstop ) and most of the energetic A = 3 elements originate from '^He reacting 
with in-flight tt^. This complicates the treatment of finding the final ^Li abundance as 
the branching to ^ = 3 elements become energy dependent. Moreover, non-absorptive 
inelastic pion-helium reactions result in mass-3 injection spectra which are continuously 
distributed. We account for all these effects with details presented in Appendix A. 3. 

In a final step we also account for potential effects on the light element abundances 
coming from the "visible" energy injection in the decays of the pions. Any primary elec- 
tromagnetic energy deposition Ei,^j is quickly dispersed in an electromagnetic cascade 
with Einj shared among a large number of photons. Photons with energies less than 
Ec — ml/22T [42] loose their ability to pair create in scatterings on the background 
radiation. Those associated "break-out" photons can then destroy the light elements. The 
earliest time at which this happens can be determined by equating Eq to the nuclear 
binding energies Ef, against photodissociation: 

_ J 2 X 10^ sec for ^Be + 7 ^ ^He + ^Ue {E^ = 1.59 MeV) 
~ I 5 X 10^ sec for D 7 ^ n p {Eb = 2.22 MeV) 

We see that the electromagnetic energy injection plays no role in the most interesting X- 
lifetime region tx ^ 10'* sec in which the lithium depleting chain (1.1) is operative. Only 
for Tx > 10* sec do we expect an influence on the element abundances, starting with the 
destruction of ^Be followed by D. The deepest bound element, ^He is only destroyed for t > 
10^ sec. Charged pions decay into /x^ which are, however, not initiating an electromagnetic 
cascade since for t > 10^ sec their rate for Thomson scattering becomes smaller than the 
muon decay rate. Thus, for simplicity, we assume that a total of E'inj = {Ee) = (7/20)m^ 
per TT^ is injected in form of electromagnetic energy, where {E^) is the average electron 
energy in the muon decay. Fortunately, though this neglects a certain fraction of accessible 
kinetic energy of the muon, a more accurate treatment is inconsequential for our further 
discussion. On the same footing, we also neglect a direct energy production in X-decay 
via e.g. X — > tt^tt^. More details can be found in Appendix A. 4. 

The pion-enriched BBN was run for different input values of ^^^^Y^, tx and the pion 
injection energy Eq. In Figure 4 wc plot an example of one possible choice with = 8, 
= 1, Tx = 10^ sec and negligible kinetic energy of injected pions. As expected, one 
can see a noticeable increase of neutrons at lithium-relevant temperatures. This leads to a 
decrease in the ^Be abundance and also an increase in D/H. The secondary source of ^Li 
although leading to a modest increase of ^Li over the SBBN prediction, turns out to be far 
below the observable level of ^Li/H ~ 10~^^. This is because the ^Li-burning reactions are 
very fast above Tg > 0.1. 

The exploration of the full parameter space yields the region in which the primordial 
lithium over-production problem is solved. For this solution, we require ^Li/H to stay in 
the interval 

■^Li/H = (1 ^ 2.5) X 10"^°. (2.37) 

The ballpark of observations lie in the range ~ 1.5 2.0 (see e.g. [4, 2] and references 
therein), but also values on the upper end of the adopted range have been reported [43, 44]. 
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^ X + n- Yx = 8, rjf = 10^ sec 
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Figure 4: Temperature evolution of light nuclei, meta-stable parent X particles and daughter tt^ 
mesons for Yx = 8 and tx = 10"^ sec as input values and for injection not too far from threshold. 
The pion induced elevation of n/p around Tg ^ 0.5 leads to the expected increase in the D/H 
and ^Li/H abundances, and to a decrease in ^Be/H when compared with the respective SBBN 
predictions (dotted lines); N is the neutron abundance normalized to baryons. 

For completeness, we also mention that the latest observations seem to suggest what could 
be called a drop-off of stars from the most metal-poor end of the Spite plateau [45, 46] 
(towards lower values of ^Li/H.) It is important to note that such a feature, as puzzling 
as it may be, in itself does not alleviate the tension between the SBBN prediction and 
the lithium observations, but rather enhances the controversy. The observational status 
on that issue is still very much under dispute, as is best illustrated by recent Ref. [47]. 
That work does not seem to confirm such a "sagging tail" in ^Li/H but rather finds two 
plateau values, ~ 1.5 and 1.9 x 10"^*^. At this point, we choose the range (2.37) based 
on the overall range of lithium abundance in the metal-poor stars along the Spite plateau. 
Should the mechanisms for astrophysical depletion of the lithium abundance firm up, the 
range (2.37) must be shifted upward. 

In addition to (2.37), we shall also put the following constraints on the remaining light 
elements which are produced in observable quantities in BBN: 



D/H < 4 X 10~^ (2.38) 

0.24 <Yp< 0.26, (2.39) 

^Li/H < 6 X 10-l^ (2.40) 

^He/D < 1. (2.41) 




-15- 



mx — 2m^± 



10^ X (D/H) 
10" X C^Li/H) 



1 



lO'' 



102 



10^ 



10^ 



Tx /sec 



Figure 5: Light element abundance yields in the plane of X-lifetime vs. X-abundance (prior to 
decay), {tx,Qt Yx) with = I and negligible kinetic energy of injected pions. The band shows 
the region 10^^° < ^Li/H < 2.5 x 10^^° in which the BBN lithium prediction is reconciled with 
observationally inferred primordial values. In the dark (blue) shaded part all limits (2.38)-(2.41) 
on the remaining light elements are respected. Contours of constant helium mass fraction Yp are 
shown by solid lines and of constant D/H abundance by dotted lines; the dashed line in the upper 
right corner corresponds to ^Li/H = 6 x 10~^^ with smaller values anywhere below. The star is the 
point in parameter space for which Fig. 4 was obtained. 

The upper limit (2.38) corresponds to the highest rehable determination [48, 49] of D/H in 
a QSO absorption system with a simple enough velocity structure. Given the significant 
scatter in the various determinations of its primordial value, the possibility of deuterium 
astration remains and values as high as (2.38) cannot be convincingly excluded at the 
present stage. For ^He the inference of its primordial mass fraction is potentially plagued 
by systematic uncertainties and values in the range (2.39) have been witnessed over the 
years with most recent determinations being on the higher side [50, 51]. Observations of 
the isotopic ratio ^Li/^Li in metal poor halo stars are extremely difficult. Though a 5% 
plateau value has been claimed [52], its inference has also been challenged [53] so that 
we only impose (2.40) as an upper limit. Finally, with D being more fragile than ^He, 
the ratio ^He/D is a monotonically increasing function of time so that the solar-system 
value (2.41) [54] provides an upper limit on the primordial value. 

In Figure 5 we present a parameter scan in the (rx, '^F^''^)-plane with sj;^^ = 1, in the 
assumption of negligible kinetic energy or injected pions, so that all reactions and decays 
occur at rest. The star is the point in parameter space for which Fig. 4 was obtained. 
The shaded band shows the region in which ^Li/H is within the observationally favored 
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Figure 6: Same as in Figure 5, but with a kinetic energy of injected pions i?o = 180 MeV which 
approximately corresponds to the maximum of the pion-nucleon cross section. 

range (2.37). More pions become available for larger values of Yx so that (^Li+^Be)/H 
decreases from its SBBN value ~ 5 x 10^^'^ at Yx = 1 to 2.5 x 10^^'^ at the lower border of 
the shaded band. As expected, there is also an associated increase of D/H from its SBBN 
value, ~ 2.6 x 10~^, with contours of constant deuterium shown by the dotted lines and 
running almost in parallel to constant lithium. The lowest D/H value along the '^Li/H 
band is 3.25 x 10~^ (for tx < 10^ sec). The contours of constant ^He are shown by the 
solid lines. The effect on ^He becomes stronger for smaller X-lifetime values. As is well 
known, ^He exhibits the strongest sensitivity on neutron/proton interconversions which 
take place before the opening of the deuterium "bottleneck" . We further observe that the 
secondary production of ^Li (2.32) becomes only important for large lifetimes tx ^ 10^ sec 
(dashed line). This, however, is already the region in the parameter space in which the 
proposed mechanism for lithium depletion by extra neutrons becomes inefficient. By the 
same token, the photo-destruction of light elements remains decoupled from the main effect 
as it only happens at late times (2.37). The destruction of ^Be and D can respectively be 
seen by the downward trend of the lithium band as well as by the dotted curve labeled 
"2". Finally, we note that the ratio (2.41) of ^He/D is never saturated throughout the 
whole (rx, ^x)-plane. The dark (blue) shading of the band corresponds to the region in 
which all constraints (2.37) and (2.38)-(2.41) are obeyed. It is this region, for which the 
pion injection provides a satisfactory solution to the lithium overabundance problem. 

As was argued earlier, pions with kinetic energy close to the delta-resonance may have 
a considerably larger efficiency of p — )• n conversion below Tg ~ 0.4. To demonstrate 
this point, we take the same assumptions as in Fig. 5 but assume a primary pion injection 
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energy of Eq = 180 MeV. The results of the new scan are presented in Fig. 6. The increased 
efficiency ofp ^ n conversion (along with stronger 7r^-^He reactions) leads to a lowering of 
the bands, i.e. less pions per baryon are necessary for the required depletion of ^Li. Note, 
however, that the overall effect is somewhat milder than initially expected. For a lifetime 
Tx = 10^ sec at which pions are still being stopped, one finds a factor ~ 2 ^ 3 reduction in 
Yx- For larger lifetimes it only becomes about half an order of magnitude (note the slightly 
different scales on the y-axcs of Figs. 5 and 6.) Though ^Be is indeed efficiently destroyed 
by strongly elevated neutron concentrations it turns out that the limiting factor is the 
final step in the overall ^Be+^Li depleting chain (1.1), i.e. the subsequent ^Li-destruction 
via thermal proton burning, ^Li + p — > ^He + ^He. This latter reaction effectively shuts 
off for Tc) < 0.3 and the overall lithium abundance is produced in form of ^Li. Residual 
late X-decays at temperatures Tg < 0.5 also lead to the narrowing of the lithium band at 
small lifetimes tx < 10^ sec and large Yx- We remark in passing that in that region, the 
accompanying increase in D also leads to an enhancement in the n-sourcing thermonuclear 
DD and DT fusion reactions. 



2.4 Decays to kaons 

Aiming at a reasonably accurate treatment of kaons is even more difficult. Let us first focus 
on K"^, and assume for simplicity that the kinetic energy is relatively small, so that one 
can neglect the effects of the decay in flight. In the case of kaon injection it is important to 

include their hadronic decays into charged pions, as well as the direct interaction of kaons 
with nucleons. The production of charged pions in kaon decays is relatively easy to account 
for: 



^±7r0(7r0) 22.4% 



5.6% 



7r^7r"'"7r 



which leads to a comparable population of and vr^. As has been already mentioned, 
there is no need to track muons (and the associated neutrino yield) since Pp^n ^ ^p-^n- 
For the average kinetic energy of produced charged pions we take Eq ~ 80 MeV. Even in 
the two-body decay, tt"^ have energies of only ~ 110 MeV so that they fall somewhat lower 
than the A-resonance, and their incomplete stopping at late times will not have a large 
effect on the efficiency of p ^ n conversion. 

We now turn to the calculations of the direct capture of kaons on nucleons. Whereas 
the charge exchange reaction K~ + p ^ + n has Q = —5.2 MeV and is thus not 
allowed kinematically on threshold, K~ reactions on nucleons can proceed via the "s-quark 
exchange reactions" with hyperons and pions in the final state: 

K-+p^ E±7r^, E°7r°, A7r°, (2.42) 
K- + n^ S-7r°, E°7r-, Att". (2.43) 



In order to obtain the inclusive n ^ p interconversion cross sections, one then has to take 
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into account the decays of the strange baryons in the final states 
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can be obtained 



(2.44) 



of an isospin (/) decomposition of the wave-functions of the reactants. Together with the 
knowledge of the relative phases ^ of the individual isospin amplitudes, the cross sections 
can then be inferred from the total if^-absorption cross sections of a given /, 



ikAj 



(2.45) 



where Ai = aj + ibj are complex scattering lengths; k is the cm. momentum of the 
incoming state. Taking into account the deviations from charge independence due to n, 
p, and K"^ and K'^ mass differences following [55, 56] and using ^0 = (—1-74 + z0.70) fm, 
Ai = (-0.05 + 7:0.63) fm, 0th = -52.9° as weh as a ratio 0.34 of A7r° to total 1 = 1 hyperon 
production [56] one arrives at the values for (2.42) and (2.43). In a final step, we use the 
branching ratios (2.44) to obtain the inclusive threshold cross sections for proton/neutron 
inter conversion, 

K-+p^n + X : {av)^~ ~ 32 mb, (2.46) 



/pn 

K- +n^p + X : {av)^~ ~ 13 mb, (2.47) 



where in (2.46) the Coulomb correction has been factored out. Whereas the first cross 
section is in good agreement with the value previously used in the BBN context, the 
second one is a factor of two smaller than in [7] . In this regard, we remark in passing that 
the latter reactions (2.43) are only due to / = 1 scatterings. 

In principle, the nucleons for some of the final states of (2.46) and (2.47) can have 
energies on the order of 30 MeV. While protons will be stopped before inducing any nuclear 
changes, neutrons with such energies are stopped primarily via interactions with p and ^He, 
and could split ^He nuclei in the collisions into some of its constituents. However, since 
^He is one order of magnitude less abundant than p, and since the maximum energy of n is 
close to the '^He desintegration threshold, we neglect such secondary effects. In addition, we 
further note that due to isospin invariance we do not need to consider hyperon production 
processes in K^n scatterings. 

As in the case of 7r~ + ^He, reactions of K~ on helium have not been considered in 
the BBN context. Fortunately, the measurements of Ref. [57] provide us with a detailed 
list of branching ratios for K~ absorption on "^He at rest. By accounting for the decay 
modes (2.44) we find the following particle multiplicities per K~ absorption 

CsHe - ^ 0.13, en = l-57-0.34AD, 

= 0.17 + 0.34Ad, Cp = 1-29-0.34Ad. (2.48) 
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When a final state is not resolved we have, for simplicity, assumed that a fraction of Ad is 
released in form of D and (1 — Ad) in form of nucleons, and we adopt Ad = 0.5. Thus, for 
example, in the reaction K~ + ^He — )■ A(S°) (pnn), which occurs with branching fraction 
22.5%, we assume 50% K{TP){Dn) and 50% A(S°)(pnn). Thereby, we are neglecting a 
certain fraction of mass-3 nuclei. In this regard, note that — unlike in the case of the 
threshold reaction ^He + vr" — an accurate computation of the secondary ^Li- yield is more 
difficult as the ^He/T injection energy is now continuous. However, as we have already 
seen in the previous section, ^Li production in excess of observationally constrained levels 
is an issue only in the region where (1.1) looses its efficiency. Thereby, we shall make 
the simple assumption that on average T and ^He carry one third of the liberated energy, 
{E^) ~ 30 McV and (-EsHe) ^ 50 MeV, respectively. (With our assumptions) the latter 
value is higher as '^He stems predominantly from processes with A and not S production. 
Finally, we note that the contribution of hyperfragments of ^He is ~ 2% [57] which thus 
does not pose any further complication. With the above multiplicities (2.48) one then 
obtains effective cross sections for each isotope in the final state, 

+ ^ iV + n : {av)'^' = ^w(^4;t'), (2.49) 

with n symbolizing an arbitrary pionic final state and N = '^He, T, D, n or p. The extrac- 
tion of {(Tv)^ is similar to the pion case and is given by (A. 8) in Appendix A. 2. 

In order to take into account the effects from electromagnetic energy injection in kaon 
decays we recall from the previous section that ir^ as well as jji^ are unlikely to initiate 
an electromagnetic cascade for tx ^ 10^ sec. The multiplicities in the K='=-decay to im^ 
(including muons from tt^ final states) and to tt'^ are 2± 1 and ~ 0.3, respectively. 
For the electromagnetic energy release in the muon decay we again take E'j'^^j ~ ra^jZ and 
for 'K^ we assume that, in addition to their rest mass, they carry on average approximately 
one third of the energy released in the particular decay channel, yielding El^^^ ~ 250 MeV. 
Taken together, this amounts to a "branching fraction" ^y^^^.^^ ~ 0.55 and a total 
electromagnetic energy injection of E\^^ ~ 2 x Br^±_^yism;^± per X-decay. As for simplicity 
we assume that the kaon kinetic energy is small, we do not have to worry about the increase 
of the electromagentic deposition by decaying in-flight kaons. However, we emphasize that 
even in the case of slow kaons, and because of the exponential sensitivity to rx, 0(10) 

per baryon will have strong effect on nuclear abundances once a photodissociation 
threshold is met for tx > 10^ sec. 

The kaon-induced solution to the cosmological lithium problem is plotted in Figure 7. 
In contrast to the case of primary pions at rest (see Fig. 5), a reduction of ^Li/H to observed 
values (2.37) is already possible for one injected kaon per baryon. This is evidently due 
to the fact that K~ cross sections on nucleons and helium are significantly larger than 
for thermal 7r~. Moreover, pions from kaon decays only amplify any effect on the light 
elements. Notice, however, that for too high of an n-abundance (large Yx) ^Li + ^Be 
starts to increase again. Though ^Be is still being depleted, net production of lithium 
in form of ''Li eventually takes over. Similar to the case of primary tt^ wc find that 
solutions corresponding to Tx ^ 10^ sec are not viable because of ^Li overproduction. 
This again excludes the potentially interesting region with tx = few x 10^ sec in which 
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Figure 7: Same as in Fig. 5, but for the case of if injection and assuming one charged kaon 
pair per X-decay, — 1, being injected close to the kinematic threshold. Because of the larger 
cross sections of on nucleons and helium, one requires a fewer number of initial X decays. The 
trend of depleting ^Be + ^Li for larger values of Yx is eventually reverted because of the enhanced 
production of lithium in form of ^Li due to numerous extra neutrons. 

previously formed deuterium is destroyed in photo-dissociation and brought back to its 
SBBN prediction. 

We now turn to the case of injected neutral kaons , — propagating as Kl and Ks- 
Whereas Ks definitely decay before interacting with the nuclear background, the relatively 
large Kl lifetime tK]^ = 5.1 x 10"^ sec may render its impact on the BBN output even 
more drastic then in the K^ case. In contrast to K^, however, Kl are not stopped by 
electromagnetic interactions so that one should take into account the energy dependence 
of their cross sections with nucleons. Again, we avoid this complication by assuming the 
injection of kaons close to the threshold. 

For i^i-nucleon scattering we exploit charge indepence and use isospin relations of the 
strong interaction processes. Indeed, conservation of strangeness implies that the inelastic 
scattering of Kl on nucleons essentially resembles the one of K^, as Kl = 2^^^^{K^ — K'^) 
up to small CP-violating corrections. Thus, up to corrections due to K^, K^, and n, p 
mass differences and other isospin- violating effects, one finds 



where Y, Y'{-ir, vr') stands for a hyperon (pion) and Y, Y'{7r, vf') are the associated states with 




a{K n — ^ Ytt) 



a{K-p Y'n') 



(2.50) 



(2.51) 
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flipped third isospin component — the final states of (2.43) and (2.42), respectively. The 
energy dependence of the cross sections not far from the threshold is inferred from (2.45). 
Indeed, since (2.50) is a pure 1 = 1 process, a{KLp — > Ytt) = ai/2. With the branching 
fractions (2.44) the inclusive cross sections for p ^ n as well as for n ^ p conversion are 
then readily found from (2.50) and (2.51), respectively. 

In addition to the reactions listed above, the following processes are possible: K^n — )• 
K~p, Klp —7- K~^n and Klp — t- Ksp- The cross sections for the former two processes are 
approximately equal, with Kj^ scattering being mediated by its or component, so 
that they can be inferred from an expression for K~p — > K^n in [56] by detailed balancing. 
Though the above processes potentially contribute to the loss of Kl and/or generation 
of K^, the probability of such processes within a kaon lifetime is very small so that they 
can be safely neglected. In addition to the processes discussed above, one should again 
account for the induced charged pion population from to and Ks decays. The respective 
branching fractions read 

Kl n+X, TT-X : 46.3%, Ks tt+tt" : 69.2%, (2.52) 

where X is the inclusive particle final state without charged pions (not to be confused 
with X-rclics). As in the case of charged kaons, the kinetic energy of outgoing pions can 
be as large as ttik/'^ — JTItt — HO MeV, and already in the regime where p + tt~ cross 
sections is enhanced. Finally, we would like to make a comment on '^Ile-K^ reactions. 
As we have just seen, Kl scattering on nucleons can be related to iC='=-reactions using 
charge independence and isospin symmetry. On these grounds there is no obvious reason 
to suspect that i^^L-reactions on helium will be drastically different in magnitude when 
compared to K~ . 

Instead of attempting to infer the ^}1c-Kl cross sections wc shall rather outline some 
general features connected to the Kl case. By turning oS the K^ reactions on helium 
in the Boltzmann code we can directly assess the difference between the charged and 
neutral kaon case. Whenever the Kl energy in the final state is very small, i.e. on the 
order of a few MeV, the results are essentially identical to Fig. 7. To conclude the kaon 
section, we would like to comment that several avenues for improvement still exist after our 
analysis. For example, the treatment of pions from kaon decays is done via assigning them a 
simplistic average energy, while of course a continuous energy spectrum would be required. 
Furthermore, the effects of incomplete stopping of K" and energetic Kl would also need 
to be included if one intends to explore the BBN sensitivity to a wider range of mx- 
Nevertheless, we also remind the reader that the effect of in-flight pions is not as drastic 
as naively expected so that the corrections due to continuously distributed vr^ from kaon 
decays will be only very moderate. Moreover, we close this section by remarking that, 
in principle, no further conceptual complications arise when considering more energetic 
injections of K"^ and Kl and that such an analysis can readily be performed by employing 
a treatment along the lines of the previous section. 

2.5 Decays to muons and neutrinos 

Stopped muons in the final state of the X-decay are sourcing energetic muon- and electron 
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neutrinos. In order to include these neutrinos in the set of Boltzmann equations we shah 
use the fact that Ts^op 

is slower than the Hubble rate. We account for the continuous 
neutrino injection and energy redshifting by calculating the neutrino phase space distribu- 
tion function f{T, En). Previous analysis of BBN modified in the presence of the energetic 
neutrinos [58] concentrated on the effects of energy deposition and neglected the direct 
nuclear-chemical impact of neutrinos [59], which turns out to be a more important effect 
at early times. 

The calculation of f{T,En) depends on the injection rate Finj, and on the primary 
spectrum of neutrinos at the time of injection, 

z/e, Ve : l2E'^{Eo - E)Eq^ for 0<E < Eq, 
Fl^{E, Eo) = \ Uf,, : 2^2(3^^ _ 2E)Eq^ for < £; < £^0, (2.53) 

for > Eq. 

where Eq ~ is the neutrino end-point energy in the muon decay, and F^^^ is normal- 

ized to unity, J F^^^dE = 1. Once injected, neutrinos are subject to flavor oscillations so 
that the primary energy spectrum will be "distorted" . Apart from vacuum oscillations, the 
neutrino propagation is affected by the coherent neutral- and charged-current interactions 
with particles in primordial plasma. To assess the importance of neutrino-refraction we 
can compare the vacuum contributions in the neutrino Hamiltonian, 

^ > 10-^^ eV, > 10-^^ eV, (2.54) 

with the matter-induced potential Vm; E < Eq. Here, Am^^j ~ 7.7 x 10~^ eV^ and 
l^"^atml — 2.4 X 10~^ eV^ are the respective mass-squared differences responsible for solar 
and atmospheric neutrino mixing [33]. For Vm one can write (see e.g. [60]) 

= ± (8 X 10-^^ eV) r| X r ^ ' ^ ' , 2.55 

[(-l + 2y,,+ye) i0YV^^r{V^.,r) 

where the Yi denote the particle-antiparticle asymmetries normalized to baryons. The 
numerical factor in front of (2.55) is Girn^/v^ with the overall +(— ) sign for viV). Unlike 
the case of charged leptons where Yf, = 0(1), the asymmetry in the neutrino sector could 
be large. In this paper we do not consider such scenarios and limit the asymmetry in the 
neutrino sector by requiring \Yi,\ < 10^. The last condition corresponds to = \fii,\/T < 
5 X 10~^ with /ijy being the neutrino chemical potential. Comparing (2.54) with (2.55) then 
implies that the flavor-evolution of injected neutrinos is given by their vacuum oscillations. 
Note that the restriction on the neutrino chemical potentials also renders neutrino-neutrino 
self-interactions unimportant which usually have the effect of locking the neutrino modes 
to each other, leading to a coherent oscillatory behavior [61]. 

A major simplification in calculation of f{T,Ei,) occurs due to the large rates for 
vacuum oscillations, 

r. osc = 1.2 X 10« C-^) (\^) > I''' ' = (2.56) 

V E A 1 eVV 1 5.5 X 10^ sec-i z = atm ^ ' 
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compared to Hubble rate and neutrino injection rate, Fj osc ^ H,Tx- This allows us to 
replace the primary injection spectrum in (2.53) by an effective one, 

Fe = {Pee)F^ + {P^e)F^ (2.57) 

where the (Pee) and {Pfj,e) are the 1/6(1^6 )-survival and v^ij^^) — )• 1/6(776 )-appearance proba- 
bilities, averaged over oscillations, 

(Pee) = 1 - ^ 26*12 ^ 0.57, (2.58) 

(P^e) = {Pe„) = \ Sin' 2012 cos2(^23) - 0.23. (2.59) 

The vacuum mixing angles arc given by sin^ ^12 = 0.312 and sin^ ^23 = 0.466 [33] and wc 
have assumed ^13 = 0. Note that (2.57) already accounts for the appropriate reduction of 
the electron neutrino flux due to Ve^pe) ^e{pe) ^t{Pt) disappearances.^ 

Defining f{T,E,^) in such a way that / f{T, E,y)dEi, is equal to the total number of 
energetic neutrinos per baryon, we arrive at the neutrino distribution function, 

where Finj = Fx, since we assume that the neutrino-sourcing muons result from the decays 
of the X-relics. 

Electron antineutrinos will have the largest effect on the BBN network. Integrated with 
the weak cross-section over energy, the distribution function gives the energetic neutrino- 
induced rate of p — )■ n conversion: 

f-Eo 

Ve+p^n + e+: V;^ = n^iT) / a-^{E,)fe{T, E,)dE, (2.61) 

where rHi{T) is the number density of baryons at temperature T. The cross section for 
quasi-elastic neutrino nucleon scattering reads (see e.g. [62]) 



a^{E,) = 0.0952 x lO'^^ j ) S{E,) cm^, Q ~ -1.8 MeV, (2.62) 



where Ei, and E^ = Ey — Aninp are the energies of Ve and in the rest frame of the 
proton, respectively; Am„p ~ 1.293 MeV is the neutron-proton mass difference and pe is 
the positron momentum. We introduce a correction factor S{Ey) = (1 — 0.0063i?i//MeV) 
which improves the agreement between the simple formula (2.62) and a precise evaluation 
of (7p„ in [63] to better than 1% in the E^-ne^me of interest. 

In our code we also account for exoergic n ^ p conversions via Vg + n ^ p + e~ . This 
process, however, is only important for tx ^ 180 sec, i.e. prior to n consumption by ^He. 
The associated cross section a^piE^) is obtained from (2.62) with Ee = Ey + Arunp and 
S{Ey) = 1. As in (2.61) we average (T^p(i?,/) over /e in order to obtain the rate FJ^^ for 
n p conversion. 



^In a similar fashion, = (P^^) F^^ + iP,^) F° and F^ = (Per) F° + {P^r) F° with / dE{Fe+F^+Fr) = 2. 
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Figure 8: Same as in Fig. 5, but for X decaying into /i^/i^ pairs which eventually decay into 
neutrinos (and e^). Note the smaller values tx in comparison to the previous figures. Energetic 
neutrinos are accumulating over time so that X-decays as early in time as i = 10 sec can affect 
the BBN network. Noteworthy is also the increased '*He abundance for large values of Yx over the 
whole rjc-range since the X-matter density is contributing to H. In this figure, m,x — 2mp. A 
novel feature is also the appearance of a second observationally favored region for tx — 10* sec. No 
additional ^Li in excess of SBBN values is produced and electromagnetic energy injection depletes 
previously formed D/H. For the parameters marked by the star the temporal evolution of the light 
elements is shown in Figure 9. 

Muon neutrions and 17^ which are sourced from muon decays (at rest) are not 
capable of interconverting protons and neutrons as their maximum injection energy 
lies below the reaction threshold 0(100 MeV). However, 1^^(17^) along with fe,r(j^e,T) are 
in principle capable of dissociating ^He via their neutral current interactions. Among the 
possible final states, only the ones with mass-3 elements have an appreciable cross section 
in the low energy regime < [64]. We include the following neutral current (NC) 

and charged current (CC) reactions into our Boltzmann network 



For the NC reactions of muon and tau neutrinos we infer the reaction rates averaged over 
ffi,T by following the same logic from (2.57) through (2.60). For the cross sections we use 
fitting formulas interpolating the tables provided in [64]. 

Given that the efficiency for neutron-to-proton conversion per each injected neutrino 
is so low [cf. (2.24)], a high initial AT-abundance (in comparison to baryons) is needed 
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Figure 9: Temperature evolution of liglit nuclei, meta-stable parent X particles and daughter ^'^ 
for Yx = 5 X 10** prior to decay and tx = 8.5 x 10"^ sec. Starting off with the previously observed 
effects which are induced by an elevated neutron abundance (suppressed ^Be/H, elevated D/H), 
electromagnetic energy injection eventually dissolves deuterium and beryllium for Tg < 0.06. This 
not only reinstates the SBBN D/H prediction (corresponding dotted line) but also implies that 
lithium is mainly produced in form of ^Li. 



in order to achieve an appreciable reduction of ^Li + ^Be. In our code we take the extra 
contribution of px to the total energy density into account. By comparing px (prior to 
decay) with the radiation energy density ftad; 

we see that a GeV-scale relic X starts contributing to the Hubble rate appreciably during 
BBN for Yx > 10^ 

In Figure 8 we show the results of our computation for X — t- p~ in the usual 

(X) (X) 

{Tx,^fi Yx) plane for =1. For simplicity, again, we assume the muon injection close 
to the kinematic threshold, mx = 2m^, although this assumption is less crucial compared 
to the hadronic case. Note that this time we have extended the tx range to lifetime values 
as small as 10 sec. This is because the usual temporal correlation between X-decay and 
its effect on the BBN yields breaks down. For example, neutrinos injected at t = 10 sec 
may still induce p — t- n conversion at t = 100 sec. As expected, the results show a similar 
pattern as in the pion and kaon cases with Yx scaled to larger values. In the neutrino 
case, however, one finds an increased ^He abundance for Yx ^ 10^ over the whole tx- 
range. This is because the X energy density is contributing to the Hubble rate, leading 
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to an earlier n/p freezeout and thereby to a higher ^He abundance. Though this makes 
results sensitive to mx , it is nevertheless only a mild dependence for nix ^ 1 GeV in the 
interesting region Yx ~ 10^, in which lithium is reconciled with observations. As it further 
turns out, the neutrino reactions dissociating ^He are giving at most a marginal correction 
to the BBN yields even for X-abundances as high as Yx = 10^. The reasons are that the 
processes of (2.63) possess threshold energies i?th ~ 20 MeV and because the associated 
cross sections are significantly smaller in comparison to (2.62). This has another interesting 
consequence: in contrast to the pion- and kaon-scenarios essentially no ^Li is produced for 

^ 10^ sec, allowing for a second region of cosmologically favored ^Li abundances in 
which also high D/H is photo-dissociated. As discussed earlier, this region is necessarily 
fine-tuned in tx because of the exponential sensitivity in injected electromagnetic energy. 
Indeed, for even larger lifetimes tx ^ 10^ sec one would observe what could be called the 
"photo-dissociation catastrophe" with ^He along with all other elements being dissociated 
(for large enough Yx)- Nevertheless, we find this scenario rather distinct if not remarkable. 
The injection of muons with relic X decays with lifetimes of 10^ seconds is capable of 
reducing ^Li abundance while keeping other abundances close to the their SBBN-predicted 
values. This occurs because of the non-monotonic evolution of D/H: first the increase due 
to neutrino-induced neutron enrichment, followed by the decrease due to energy injection. 

The temporal evolution of elemental abundances for a single point in the parameter 
space marked by the star in Figure 8 is illustrated in Figure 9. The most prominent feature 
is the "bump" in the neutron abundance at Tg ~ 0.06 marking the point in time in which 
deuterium is subjected to photo-destruction. The associated decrease in D/H from its 
elevated values down to its original SBBN prediction (dotted lines) are clearly seen in this 
Figure. Also ^Be is dissociated around the same time so that the overall lithium abundance 
comes essentially in the form of ^Li. 

3. Metastable GeV-scale states and the lithium abundance 

In this section we consider some specific models of GeV-scale relics and their impact on 
the BBN predictions. All models can be subdivided into broad categories of WIMPs 
and super- WIMPs, where "super-" refers to the superweak interaction strength with SM- 
states. WIMPs are thermally excited above their mass scale, then depleting their number 
density via annihilation at T < mwiMP, and — in our case — decaying after the start of BBN. 
Superweakly interacting particles typically have very small production rates throughout the 
whole history of the Universe. This justifies the assumption of super- WIMPs being initially 
absent as the Universe enters its thermal radiation-dominated stage e.g. after inflation and 
reheating, so that only some small super- WIMP abundance develops due to the thermal 
leakage from the SM states. Let us stress again that the notion of a (super-)WIMP here 
is not to be confused with a dark matter state (which has to be stable on cosmological 
timescales.) Of course, some of metastable particles presented in this section are motivated 
by their potential as mediators of the interactions between dark matter particles and the 
SM sector [27, 65]. 
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The simplest WIMP model without any need of UV completion [66, 67] uses the so- 
called Higgs portal. In this model, a singlet scalar field S interacts with the rest of the SM 
via its coupling to the Higgs doublet H, 

>CH-portai = lid^Sf - V{S) - {XSS + AS){H^H). (3.1) 

Here we assume an approximate Z2 symmetry, S ^ —S, broken only by the last trilinear 
term AS{H^H), and we further take A to be very small; V{S) is the scalar potential in 
the secluded sector. If the A-coupling is above O(10~^°), the S'-particles are guaranteed to 

be in thermal equilibrium with the SM-states as soon as the plasma temperature is around 
the electrowcak scale. Smaller values of A make S* to a super- WIMP [68]. The physical 
mass of S particles comes from the potential term m'^S'^ /2 and the electroweak vacuum 
expectation value (VEV) v ~ 246 GeV, m| = + \v^. Having 

A, A, and m| (S-portal) (3.2) 

at our disposal, we can always choose the region of parameter space where the lithium 

abundance is reduced along the lines described in the previous section. Of course, in this 
model the GcV scale is not special, and metastable S'-particles at the electrowcak scale can 
also be used for the same purposes, exploiting nucleons in the decay products of S. 

Going away from the simplest possibility, we introduce a U(l)5 vector portal model [69], 

>Cv-portai = -\v^^. - \ F^V^^ + \D,ct>\^ - V{<t>), (3.3) 

where the connection between the U(l)5' field strength V^y and the hypercharge field 
strength is mediated via the kinetic mixing parameter k. The U(l)5 covariant deriva- 
tive is given by Dfj_ = + ie'V^, with being the new vector-state associated with V^i, 
and e' the gauge coupling strength in the secluded sector. 

Since we are going to consider GeV-scale phenomenology, one can substitute the hyper- 
charge with the photon field strength, F^^, F^y and absorb the cosine of the Weinberg 
angle into the new definition of k. After the spontaneous breaking of the U(l)5 gauge 
group by a Higgs' field ^, the low-energy Lagrangian can be written as 

>C = -\vl + \mlV'^ + \{d,h'f - \mlh'' + Ant, (3.4) 

where my = e'v' becomes the mass of and rrihi is the mass of the physical Higgs' field 
/i', {(f)) = f'/\/2. Assuming a standard Higgs potential in the U(l)5 sector, the interaction 
terms are given by 

22 22 

Lint - "2 V-f* + V + ^ V - 2^/* - • i-^-^) 

Thus, like in the previous example, the model is characterized by only a handful of 
free parameters: 

a', K, mfii, and my (V-portal), (3-6) 
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where, as usual, a' = e'"^ /Att. The strength of the mixing angle k is undoubtedly a very 
important parameter, as it allows to make the link between the SM- and U(l)5-sector 
arbitrarily weak. Indeed, for n < 10~^^ any production rate of h' and V particles is smaller 
than the Hubble rate, and the model becomes a good candidate for the super-weak regime. 
It is also important to notice that long lifetimes of particles from the U(l)5 sector can 
be achieved without requiring exceedingly small k. This happens rather naturally in the 
regime nih' < rny, where the decay amplitude of the Higgs from the secluded sector is 
suppressed by k? [32, 70]. As mentioned in the introduction, the (sub)-GeV mass scales 
for 11(1)5 states are motivated by the possible enhancement of annihilation of weak-scale 
WIMPs into leptons [23, 24]. 

3.1 WIMP regime 

Here we would like to determine the values of parameters for both models, (3.2) and (3.6), 
that lead to the desirable depletion of the ^Li abundance. Since in the WIMP regime the 
coupling constants are not necessarily small, the results of this subsection might be relevant 
for the direct searches of new physics at GeV energy scales. 

For the Higgs portal model in the WIMP regime there are several generic choices 
of parameters that lead to the desirable range of the lifetime-mass-abundance "islands" 
that suppress the overall lithium abundance. These regions can be readily found using 
the results for the lifetime [71] and freeze-out abundance of ^-particles [67, 72, 73]. For 
example, the following choice of parameters, 

Ys 0^3x10^, rx~ 500 sec, = 250 MeV =^ ~ 2 x 10-^ A = 4 x 10"^ GcV, 

(3.7) 

leads to the depletion of lithium via the decays of S-particlcs to muon pairs, generating 
electron antineutrinos. Higher masses of S particles, ms > 2m7r that generate hadrons 
in the decay product would typically require A ~ 0(1) in order to have Ys ^ 10^. The 
estimate (3.7) is performed for a Higgs mass value of 100 GeV but can be easily rescaled for 
heavier values. Notice that the A parameter does not enter in the abundance calculation, 
and can always be adjusted to obtain a desirable lifetime. 

Interestingly enough, the mass range of ms ^ 2 GeV is being probed through the 
search of missing energy decays of S-mesons, B K(*^SS, [72, 73, 74] as the eventual 
decays of S particles occur far away from the 5-decay vertex. When the phase space 
suppression can be neglected, the prediction for the missing energy branching ratio at 
mh = 100 GeV is Br^^^^ 4 x 10"^ -|- 3 x 10"^ A^ [72], where the first term stands 
for the SM contribution via B — t- Kuv'. The current upper bound on the branching ratio 
stands at 1.5 x 10^^, excluding models with A ~ 0(1). The choice of parameters (3.7) is 
still (barely) allowed, both by the B and K decays with missing energy. We conclude that 
only this low ms option (or alternatively ms ^ 2 GeV) is capable of solving the lithium 
problem without enhancing the missing energy branching ratios of B-mesons beyond what 
is observed. Perhaps the most interesting consequence of large values for A required by the 
solution to the lithium problem are the implications for Higgs physics. For Higgs masses of 
150 GeV and lighter the missing energy decays h — >• SS will dominate over the Standard 
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Model width and lead to a significant enhancement of the missing energy channel and 
suppression of "visible" decay modes [67]. In short, the viability of the proposed scenario 
will be tested at the LHC, (super-)5-factories and new kaon facilities. 

Metastable GeV-scale particles from a U(l)5 sector is another interesting possibility. 
In the WIMP-regime we are drawn to consider h! since the vectors V decay well before 
BBN through their now appreciable kinetic mixing k, S> 10~^° [see Eq. (3.18) below], 

^<„.05.ecx(i^)'(555J^) for my>n,,. (3.8) 

Indeed, the longevity of h' particles can be achieved rather naturally in the regime m/j/ < 
my- The lifetime of the h! particles due to the four-body decays and one-loop induced 
amplitudes was calculated in Ref. [75]. For example, the decay width to muons in the 
regime 2m^ < rrihi <C my is dominated by loop effects and is given by 

, , 4, /a\ /3.4x lO-^X"^ /250 MeV\ / my V /o „x 

Th' ~ (10^ 10^) sec X — ^TT^Trv^ • (3-9) 

" ^ ^ Ka'J V K J \ rnh' J V500 MeVJ ^ ' 

This formula remains approximately valid in the regime m^' ~ my, and breaks down at 
my — m^' ^ m/j/. With my ~ 500 MeV and a' ~ a it is easy to see that for m^/ = 250 
MeV and k ~ (3 ^ 10) x 10~^ the decay rate of ^' is in the right ballpark, while for the 
same value of couplings and masses the decay of vectors happens very fast, ry < 10~^^ sec. 

Since we take the h! particle to be the lightest in the U(l)5 sector and long-lived, 
what physical mechanisms could possibly deplete its abundance to an acceptable level? It 
turns out that possible endothermic excitations into the V state, h! ^ V followed by V 
decays, although exponentially sensitive to the mass parameters in the U(l)5 sector, are in 
principle sufficient for the h! depletion^. In particular, the following processes are capable 



of suppressing the h! abundance: 

h' + h' + V, Ti oc (a')^K°exp(-m^//r-2Am/r) (3.10) 

h' + V^l-^r, r2(xaWexp(-m?,//r-Am/r) (3.11) 

h' + 1^ + Ts (X a'aK^ exp(-Am/r), (3.12) 



where Am = my — m/j/, and / are the charged particles of the SM. The last process 
in (3.12) is especially important because it comes with the least amount of exponential 
suppression. We now estimate the freeze-out abundance of h' due to this process. 

In the norelativistic regime, the equilbrium number density n^^ of h' follows the well- 
known curve, and the freeze-out abundance can be estimated by equating the depletion 
rate with the Hubble rate, H{Tf) = T^iTf), 

Yh,^n'%Tf)/n,{Tf), (3.13) 

= (^) «^p(-^'^'/r). (3.14) 

^MP would like to thank Neal Weiner for a very useful discussion of that point. 
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We would like to determine the dependence of the freeze-out temperature Tj on Am, 
and specialize it to two cases considered in the previous section: decays to muons and 
decays to pions/kaons. Considering the exponential dependence on the mass splitting Am 
we are making a series of simplifying approximations: mg <C <C Am <C my,m/j'. In 
that case the cross section for the h' + e"^ ^ V + e"^ process is given by 

a = 4 ^ =^ rg T ~ ^ ^exp -Am/r , 3.15 

rriy he TTUiy 

where Eg is the electron energy. For Tf larger than 100 MeV, other charged SM states 
would have to be included in the coannihilation process (3.12). However, since we are in 
need of a suppression of the h' abundance by a large factor, the freeze out of a GeV-scale 
m'^ has to happen below 100 MeV, so that only 6=*= will contribute to depletion process. 

Introducing two dimensionless variables Xf = rrih' /Tf and 6 = Am/m/^/ = mv/mh' — 1, 
we arrive at the following analytic estimate of the freeze-out temperature and the resulting 
abundance as a function of mass splitting in the U(l)5 sector: 



X/ ~ - < 13 + In 



a 



(3.16) 



a (1 + 5)4 Xf \3x lO-sy nih' 
Yh' ~ 4.2 X 10** X xf^expi-xf). (3.17) 

It is easy to see that for mass splittings (5 2± 0.3 + 0.5 the abundance varies in the interval 
10~^ + 10^, which will be suitable for the solution of the lithium problem using the decays 
of h' to pions and kaons. Mass splittings of order S c± 0.8 -v- 1.4 correspond to abundances 
on the order 10^ -i- 10^, which are suitable for suppressing lithium via the decays of h' to 
muons. We remark in passing that (3.17) represents the fraction of h' relative baryons after 
-annihilation but prior to their decay. 
Since the realization that a new attractive U(l)5 gauge force can be used for the 
explanation of the PAMELA anomaly [23, 24], a lot of dedicated work has been done in 
order to understand the prospects of searching/detecting particles from a putative U(l)5 
sector. The most promising avenues for the discovery of new GeV-scalc particles arc high- 
luminosity medium-to-low energy experiments [75, 76, 77], although some prospects of 
discovering such particles in the fragmentation of heavy exotic particles produced in the 
high-energy collisions have also been investigated [78]. The mixing angles deduced in this 
section, k ~ 3 x 10~^, and vector masses in the range of 400 MeV and larger, represents 
one of the most challenging corners in the my-K parameter space. These angles are small 
enough so that the search of V in this range at even the highest luminosity e~^e~ machines 
is not possible. Therefore, the search for V in this regime with fixed target experiments is 
perhaps the only realistic option. Vector particles are relatively short lived, cr < 1 cm, so 
that in the set-up with a detector at some macroscopic distance behind the target all the 
decays will happen before reaching the detector. With proton beams, where the detector 
is typically tens or hundreds of meters behind the target, such a search might prove to 
be very challenging. Therefore, the best discovery potential for V in this parameter range 
would probably be a high- intensity electron beam on a thin target [79] . Lastly, the lifetimes 
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of mediators in excess of milliseconds lead to interesting effects in the annihilation of dark 
matter, when the decays of mediators occur away from the point where the annihilation 
occur [80, 81, 82]. This leads to novel signatures in the indirect detection of dark matter in 
models with light mediators. Finally, for over density constraints on some variants of these 
models see, e.g [83]. 

3.2 GeV-scale super- WIMPs 

In the previous examples we considered metastable particles that initially had thermal 
abundances. This is not the only possibility, and in the following we will look at the 
abundances of vectors V and scalars S in the super- WIMP regime. Our assumption is 
that the link between the SM and the secluded sectors is extremely weak, and given by a 
single parameter, n for (3.3), and A for (3.1). The states X = S or V are then produced 
with sub-Hubble rates from the thermal scattering of SM particles. In this case, both, the 
production and the decay rates are proportional to the square of the small parameters n 
or A. Consequently, the product of abundance and lifetime, Yxtx is independent on A 
or K and is only a function of the mass mx and the SM couplings/masses. This scaling 
holds for both cases despite the fact that the physical prodiiction mechanisms of S and V 
are vastly different. An extra vector particle that appears just as another massive photon 
is mostly produced at temperatures comparable to its mass T > my whereas a scalar 
particle with mixing to the Higgs boson is most efficiently produced at T ~ niw, where 
the T/m/i-dependence of the production rate relative to the Hubble rate is maximized. 

We begin by analyzing the U(l)5 model (3.3) in the super- WIMP regime. The decay 
widths of the vector particles to leptons and hadrons are given by [75] 



where R is the experimentally measured ratio of the inclusive hadroproduction in e~^e~ 
collisions to the direct muon production, R = (7e+e-->.hadrons/'^e+e-->^t+^j- • It follows im- 
mediately [see also Eq. (3.8)] that a desirable lifetime range ry > 0(100 sec) requires k in 
the ballpark of IQ-^^ 

The production of GeV-scale F-bosons in the early Universe cannot be calculated 
"exactly" because of the unsurmountable difficulty in treating hadrons in the intermediate 
regime T ~ Aqcd- In particular, all production modes involving free quark and gluons, 
qq^V,q + g^q + V etc are not calculable because of the strong QCD coupling. 
Cases that can be reliably calculated are my <C 100 MeV, and my 3> 1 GeV, which are, 
unfortunately, exactly opposite to the regime we are interested in. Nevertheless, we can 
obtain an order-of-magnitude estimate that captures the main behavior of the production 
mechanism as a function of my. The main feature of the production mechanism is that 
it receives an exponential cutoff exp(— my/T) in regime T < my, and this property is 
not modified by strong dynamics. Therefore, the production effectively stops at T ~ my 




R{s = my). 



(3.18) 



(3.19) 



-32- 



with the residual abundance of V particles parametrically dependent on the ratio of the 
production rate over the Hubble rate, 



rib s 



. (l..,4.9, X ^T-^'ir^ - 0.3 . i'-!^) (9^r (^f 

nbKs{mv) H{mv) \ ry J \mv J \gesj 

(3.20) 



The estimated range follows from a computation of the collision integral for the inverse 
decays of Fermi-Dirac distributed pairs (lower value) and massless e^, /j,^, uu, dd, ss 
(upper value); for a related calculation see also [84]. The upper value is saturated for 
"freeze-out" temperatures above the QCD hadronization scale, Aqcd < Tf < 0{1 GeV), 
where the perturbative treatment of the light quarks becomes accessible. On the other 
hand, lower mass vectors whose production ceases later are guaranteed to receive a con- 
tribution to their abundance from electron-positron pairs. This limits Yy from below. In 
the last relation we have normalized the effective number of entropy degrees of freedom 
heS ^ geS on a typical value during the QCD hadronization epoch and further made the 
assumption that ry is dominated by the channels which are also responsible for the inverse 
decay. 

The results suggest that a sub-GeV state with a lifetime of the order of a thousand 
seconds is abundant enough to produce 0(1) pions/kaons per baryon, and thus capable 
of reducing the lithium abundance by a factor of a few. In the most optimal range with 
iTiv ~ rup, the pion branching dominates over lepton branching, while the energy of injected 
pions is close to the delta-resonance. This enhances the the efficiency of p ^ n conversion. 
Fig. 6, so that Yy ~ 2 is required. Since Yy is not calculable exactly, it is unfortunately not 
possible to give a final judgement whether the required abundance can be achieved or falls 
marginally short. We also notice that my below the di-pion threshold, 2m^ < my < 2m7r 
can only be created in abundances less than 0(10^) which renders this parameter range 
incapable of reducing lithium via F-decays into muons. On the other hand, a heavier mass 
range for a super- WIMP V is also interesting. For example, F-bosons in the mass range 
of 10 GeV or more are less abundant than protons, but will contain nucleons among their 
decay products. In this case one can easily suppress the lithium abundance via direct 
injection of extra nucleons along the lines of Refs. [7, 8]. 

The model with the singlet S mixed to the Higgs boson via the ASH^H coupling can 
be treated very similarly. Interestingly, it turns out that in this model the abundance of 
S particles can be calculated without significant QCD uncertainties because it dominatly 
occurs at electroweak scale temperatures, whereas the lifetime of the S boson with a mass 
in the GeV range is notoriously difficult to handle. 

In principle many SM scattering processes may contribute to the emission of the 5 
boson. These include /1/2 ^ VS, fiV ^ /s^, hh SV, VV SV, tg tS, tt gS 
etc, where /i,/2 stand for the SM fermions, V = W,Z for massive gauge bosons, and g, 
t are gluons and top quark. An exact treatment of the production mechanisms in this 
model falls outside the scope of this paper. However, despite this rather large number 
of production channels, the asymptotic behavior of the production mechanisms are quite 
obvious. In the regime T ^ v the production rate has to scale as F oc A'^/T, and in the low- 
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energy regime it is oc A^T^rri^'^. This implies that V is naturally peaked at temperatures 
around the electroweak scale. 

In order to avoid straightforward but tedious calculations of the S abundance in the 
most minimal model (3.1), we go to the two-Higgs doublet model analogue of (3.1), and 
couple the S scalar to a "mixed" Higgs portal: 

ASH^H ^AS{HiH2 + h.c). (3.21) 

Assuming a mild hierarchy of scales in the Higgs sector together with a ratio of the Higgs 
VEVs tan^ = V2/v\ S> 1, i?2 represents the SM Higgs doublet while Hi contains heavi(er) 

physical scalars if, A, with common mass scale mn- The production of scalars S 
is peaked around mu- In what follows we calculate the abundance of S resulting from 
Hi — )• H2S decays using mn 3> mh,miY(^z)- exact analogy to (2.10) the decay rate can 
be written as TH{mH /E), where 

is the decay rate of Hi particles at rest. Accounting for that Hi carries four degrees of 
freedom we arrive at the following integral formula for the freeze-out (or rather " freeze-in" ) 
abundance of S bosons realtive to baryons. 



Ys = ^ 
rib 



J H{T)TsiT) J E(27r)3 exp(£;/r) - 1 ^ ' ^ 



To good approximation we take hes{T) in s = /4:5 to be that of the SM, heS = 

427/4. The integrals in (3.23) can be taken exactly, resulting in 

= il^-if^. 3.8 X 10= (3.24) 

where H{mH) is the value of the Hubble constant at T = mn and C,{x) is the Riemann 
Zeta-function. 

The decay of the S-bosons at late times proceed to muons and hadrons depending on 
what channels are kinematically allowed. At ms <^ 1 GeV the decay rates can be obtained 
within the framework of chiral perturbation theory and by low-energy theorems [85, 86]. For 

^ 1-5 GeV or so, pcrturbativc QCD starts getting applicable. Around these energies 
the decay rate is dominated by S decaying into s-quarks with a ~ 25% contribution to 
muons: 

Notice that heavy-quark-mediated S ^ gg decays are suppressed relative to the SM by a 
factor of 9 because in the large tan (3 regime only 6-quarks contribute, and t, c contributions 
are subdominant. Since rris ^ m^, the decays to strange-antistrange pairs happen in at 
least 75% of the cases. When ms is below the di-hyperon threshold, each of these decays 
lead to either K~ or Kj^ in the final state. The last step in checking the suitability of this 
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model involves trading in the formula for the abundance (3.25) for the lifetime of S. 
Taking nis ~ 100 MeV, we get 

where all parameters are normalized on their "natural" values. It is in a way remarkable 
that in these models many numbers of very different orders of magnitude conspire to yield 
abundances in the right ballpark necessary for a suppresion of the ^Li abundance via the 
emission of kaons (and pions) . 



4. Conclusions 



The current discrepancy between the prediction of the lithium abundance in SBBN and 
its observation at the Spite-plateau value prompted a number of particle-physics related 
explanations. Previous works have mainly concentrated on the effects of heavy, electroweak 
scale particles on the ^Li abundance. Such models typically involve decays of electroweak 
scale WIMPs with nucleons among the decay products, or the catalytic suppression of '''Li 
by metastable charged particles, that have their masses again in the electroweak range. 
In this paper we have considered a related but nevertheless distinct possibility: decays of 
(sub-)GeV-scale neutral relic particles. Because of their kinematic constraints, neither of 
these particles can be charged nor can their decay products contain nucleons in the final 
state. Nevertheless, the reduction of the lithium abundance may occur simply because 
pions and kaons in the final state have a significant probability of interacting with pro- 
tons and creating a neutron-excess. We have shown that 0(1) — 0(100) abundances of 
metastable particles with respect to baryons and with lifetimes around a few hundred to 
ten thousand seconds can easily reduce the lithium abundance by a factor of a few and still 
be consistent with measurements of primordial deuterium. We have also shown that if the 
decays to the charged hadrons are kinematically not possible, the muons in the final state 
can also reduce the overall lithium abundance. In this case the main effect comes from 
electron antineutrinos that induce p — t- n inter conversions. It turns out that O(lO^) muon 
decays per proton around a few hundred seconds are needed in order to reduce the lithium 
abundance to observationally favored levels. Moreover, injection of pions and muons not 
far from thresholds also leave less "damage" in form of ^Li, compared to scenarios with 
electroweak-mass decaying relics. This is because the ^Li production requires significant 
amounts of energy to split ^Hc and create energetic ^ = 3 nuclei, which in case of the 
muons/pion injection is far less efficient. This broadens the acceptable lifetime range to 
up to O(IO^) seconds. 

We have done a careful numerical investigation of this problem by including charge 
exchange and absorption reactions of pions, kaons and (anti)neutrinos on nucleons and ^He. 
In doing so we have updated the values for the most relevant reactions, and integrated them 
into a Kawano-based nuclesynthesis code. We have also accounted for the change in the 
neutrino spectrum due to the Hubble expansion and due to neutrino oscillations. The 
results are presented in the form of abundance-lifetime plots. Figures 5-8, where we have 
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assumed a maximal branching into each species. We can see that the elevated deuterium 
abundance is an inevitable consequence of all such scenarios, resulting from the "removal" 
of extra neutrons by the np — > dj reaction. The only exception is the case of relic particles 
decaying to muons at ~ 10^ seconds. In this case, the late energy injection, combined with 
early neutrino-induced p — > n interconversion can lead to the SBBN value of deuterium, 
and achieve the suppression of ^Li by a factor of a few. 

In order to relate this proposal to some concrete model realizations, we have considered 
two minimalistic extensions of the SM. We chose to extend the SM by a singlet scalar S 
coupled via the so-called Higgs portal, or a new U(l)5' group that has kinetic mixing with 
the photon field strength. Both models, in principle, possess new long-lived states, and 
in the case of the U(l)5 model this is either the vector V or the Higgs' particle h' of 
the U(l)5-breaking sector. The cosmological history in these models then depends on the 
strength of the couplings that connect the secluded sectors to the SM particles. 

If the couplings to the Standard Model are large enough, the new GeV-scale states 
X are thermally excited (WIMP regime). An acceptable abundance of X = S, h' prior 
to decay is then either achieved via annihilation, e.g. SS ff, or via the excitation to 
the unstable state, h'e — t- Ve — > eee. We have shown that in the WIMP regime both 
models provide viable solutions to the lithium problem, and that both, the vr, X-mediated 
or the antincutrino-mcdiatcd suppression of ^Li indeed work. We have identified the plau- 
sible choice for the parameters where such reduction will occur. The same parameters 
may be chosen to fit the PAMELA and FGST anomalies, when the mediator models are 
complemented with electroweak-scale WIMPs. In addition to cosmological probes, there is 
a certain chance of detecting these new GeV-scale states in laboratory experiments, such 
as in rare decays of mesons, in Higgs decays, and in fixed target experiments. Last but 
not least, we remark that a nih' < my mass pattern naturally arises in straightforward 
super symmetric generalizations of the ^^(l)^ model. 

We have further shown that the super- WIMP regime of these models also have a 
potential to provide a viable resolution to the lithium problem. In this case, the metastable 
states are never in thermal equilibrium and are only produced via an extremely small 
thermal rate. Such models are relatively predictive because the mass and the coupling 
of X can be translated directly into lifetime- and abundance-predictions. Unfortunately, 
the extreme smallness of the coupling constants prevents direct experimental searches of 
such particles in the laboratory. On the other hand, the sensitivity that the early Universe 
exhibits to new physics is once again emphasized by the fact the lithium and deuterium 
abundances are capable of probing effective coupling constants as small as aK^ ^ 10~^^. 
It remains to be seen whether the ^Li reduction can be achieved in another well-motivated 
model of super- WIMPs in the form of heavy sterile neutrinos. 

To close up, the problem of the ^Li abundance definitely deserves special attention, 
as it hints at one of the very few inconsistencies in the standard cosmological picture. Of 
course, at this point it is difficult to insist that particle physics must necessarily be key for 
solving the current lithium problem. Yet, as shown in this paper, several particle physics 
scenarios that reduce the '^Li abundance are quite natural and perhaps deserve further 
detailed considerations. As to the lithium problem, only future progress in cosmology, 
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astrophysics, particle, and nuclear physics may help to resolve this intriguing problem. 
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A. Appendix 

A.l Some details on tlie Boltzmann code 

The Boltzmann code that we use is based on the one of Ref. [87], in which we incorpo- 
rate some significant improvements and updates: physical constants, isotope masses, and 
conversion factors arc determined from the evaluations [88, 89]. For all important SBBN 
reactions (i.e. up to A=7) we employ the results of [90] with the exceptions of the n{p, 7)D 
and '^He(a,7)^Be reactions for which we follow [91] and [92], respectively. All other 0(80) 
remaining reactions we update by following the recommendations of [93]. 

In order to arrive at an accurate prediction of ^He we numerically integrate all weak 
rates at each time step for which zero temperature radiative corrections are taken into 
account following [94]. When applicable, we further take into account relativistic Coulomb 
corrections to the weak rates by multiplying the integrands with the Fermi function 

F{Z,v) = 2(1 + S)i2pR)-'^^-')\^!j^^^e-^ (Al) 

where Z is the nuclear charge, S = {1 — o?Z'^)^l'^ ^ p (v = p/E) is the momentum (velocity) 
of the relative motion, and i? ~ 1 fm is the proton radius; 77 is given in (A. 4). In our code 
we approximate F{Z, v) by [95] 



F{Z, v)^Y^ an{aZY, with 

"11 



n=0 



vr 11 1 TIN '^^ ^ 

ao = 1, ai = -, 02 = — - 7E - In {2pK) + — ^, 03 = - 

V 4 iv^ V 



^ -fE-\ni2pR) 



(A.2) 



where — 0.5772 is the Euler-Mascheroni constant. Note, however, that in the non- 
relativistic case, F{Z,v) is better approximated by (A.3). As the dividing line we choose 
V = 0.045 which guarantees that (A.2) and (A.3) represent (A.l) to better than 0.1%. In 
the non-relativistic limit the familiar factor reads 

where Z is the nuclear charge; rj is called the Sommerfeld parameter. Evidently, note 
that for Coulomb corrections of equally charged particles one has to replace Z by —Z, 
e.g. in 7r++^He reactions. When obtaining Coulomb corrected cross sections, we make the 
approximation 

{av)i ~ F{Z, {v))iav)i, {v) = ^JWf^ (A.4) 
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where ^ is the reduced mass of the incoming particles. 

Finally, as assessed in [94], we apply a slight upward-shift of 0.72% of the resulting ^He 
abundance in order to account for more subtle, subleading corrections to the BBN reaction 
network. We then find very good agreement of our light element abundance predictions 
with the ones presented in the recent work [3] at the WMAP value of rji, = 6.23 x 10"-'^° 
and with a neutron lifetime of r„ = 885.7 s. 

A. 2 Pion and kaon capture on helium 

Experimental data on absorption of stopped charged mesons on various target materials 
are usually cast in the form of absorption rates and branching fractions of the reaction 
fragments. After being stopped, ir~ are captured into atomic orbits from which they emit 
X-rays (and Auger electrons) until they come close enough to interact strongly with the 
nucleus. The cascading time is significantly smaller than r^± and pions reach a IS orbit 
prior to absorption. A measurement of the ground state level width Ti^ = (45 ±3) cV [36] 
of pionic helium then allows us to obtain the low-energy in-flight cross section (av) of a 
free ^He-7r~pair via 

^fs = \i^ism\<Tv), (A.5) 

where |V'i5(0)pv is the flux density at the origin. Assuming a 100% absorption from the 
IS state, we find for the total 7r~-^He reaction cross section 

(ctIh^u) ~ 7.3 mb. (A.6) 

To extract the cross sections for the K~ + ^He reaction we can proceed in a similar 
fashion as for vr" -|- ^He. Interestingly, however, measurements of kaonic X-rays from ^He 
show no radiative transitions from the 2P — t- IS* state [96] which tells us that on the 
2P-level prefers to react strongly with the alpha nucleus rather than undergoing another 
radiative transition to the ground state. Since ~ 16% of captured K~ indeed make it to 
the 2P-state, this implies a lower bound on the 2P absorption rate [96] 

r|^^ > 3 X 10^^ sec-^ (A.7) 

Using that value will, however, only yield the (momentum-dependent) p-wave contribution 
to the total in-flight absorption cross section. Therefore, and since we could not locate 
further experimental data on s-wave absorption, we have to rely on theoretical calculations. 

Reference [97] finds for the absorption rate from the 2S'-state r|^^ = 48.1 x 10-*^^ sec"-*^ 
which points towards a very large cross section (av) 2± 1.2 b [using 'ip2s instead of tpis 
in (A.5)]. This value, however, is in conflict with the s-wave unitarity bound (cJmax^) — 
0.4:h/\/Jg so that we refrain from using it. We further observe that the quoted value 
for the 2P-absorption rate in [97] is a factor of ~ 20 larger than the experimental lower 
bound (A.7). Employing the (somewhat ad-hoc) prescription to scale down of [97] by 
the same factor 20 (which saturates the lower bound for the 2P-case), we arrive at a total 
in-flight threshold cross section of 

(a^^^v) ~ 60 mb. (A.8) 
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Figure 10: Left: Cross section for pfon-induced proton conversion, tt~ + p ^ ir" + n, across the 
A'^-resonance as a function of pion kinetic energy in the LAB-frame. The sohd Une shows a pion- 
nucleon partial- wave model calculation [99, 100]. In-flight experimental data labeled "Brs06" are 
from [101], the point at threshold labeled "Gas08" is inferred from the IS level-width as quoted 
in [34]. Right: The top curve shows the total 7r*-^He inelastic scattering cross section as given 
in [38] and labeled as "Bin75". The data labeled "Mat98" [40] is the absorptive part; older data 
"Str9G" on pion absorption is collected in [39]. We employ the model calculation of [102], scaled 
by a factor 0.8 to fit "Bin75". The dashed line shows the difference between the latter two lines 
and thus provides the remaining part of the inelastic cross section. The dotted parts of the lines 
indicate the regions in which we extrapolate without data. The dot at threshold is again inferred 
from the pionic helium IS level- width given in [36]. 

Indeed, this cross section compares well with the values inferred from another calcula- 
tion [98] for absorption from the IS-state which gives us some confidence that (A. 8) 
lies in the right ballpark. 

A. 3 Pion reactions across the Delta-resonance 

The predominant feature in low-energy pion coUisons with nucleons and nuclei is the forma- 
tion of A-resonances. For reactions with nucleons, we only need to consider ■it~ +p — )• ir^+n 
across the resonance. Though at at cosmic times t < 180 sec neutrons are still an abundant 
target, tt^ are stopped before interacting with them. For the charge exchange reaction of 
TT~ +pwe employ the partial wave analysis of [99, 100] as shown in the left panel of Fig. 10 
as a function of pion kinetic energy E.^ in the rest frame of the proton. It agrees well 
with the experimental results of [101]. The dot at threshold corresponds to (2.26). The 
alternative exit channel n -|- 7 is only relevant at threshold and becomes quickly negligible 
with increasing pion energy. We have verified this from a model calculation [99, 100] of the 
inverse process of photo-pion production by employing the principle of detailed balance. 
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Table 1: Adopted final state multiplicities for the various 7r**He reactions. Values are rounded 
such that baryon number remains conserved, 3^3He + 3^t + 2^d + + £,p = 4:. 

The right panel of Fig. 10 shows the various contributions to 7r^-^He scattering. The 
top curve and data "Bin75" represent the total inelastic cross section as taken from [38]. 
The cross section breaks down into an absorptive part with no pion in the final state and the 
remaining inelastic part. We fit the in-flight absorption data "Mat98" [40] by scaling down 
the calculation of [102] by a factor of 0.8; older data "Str90" on absorption is collected 
in [39]. The dotted parts of the lines show the regions in which we had to extrapolate 
without data. In particular, note that inelastic pion scattering (without absorption) has 
thresholds in the range ~ (15 30) MeV depending on the final state. For simplicity, and 
because absorption dominates for low energies, we have extrapolated the inelastic cross 
section to a common threshold value E^^^ ~ 25 MeV. 

In particular, the output of ^He and T from absorption and inelastic vr^-'^He scatterings 
is important for the non-thermal ^Li abundance. The adopted nuclear /nucleon final-state 
multiplicities are summarized in Table 1. For absorption, one finds that the branchings 
into the various final states differ significantly depending on whether the process occurs 
at threshold or not. At threshold the fraction of ^ = 3 elements is ~ 20% but for E^^ = 
120 MeV the respective multiplicities, as inferred from the measurements of [39], are already 
very small. Since, by nature, absorption yields the largest recoil energies such a behavior 
is not unexpected and between E^^ = and 30 MeV we interpolate between the threshold 
and in-flight ^-values. For larger pion energies we switch to the in-flight data. 

Above the threshold of inelastic 7r^-^He scattering the energy-dependent multiplicities 
of elements i = ^He, T, D, n and p are given by 

CSel = fsxClIx + (1 - /sx) [(1 - 9DX)^I^X + 9DX^^.Dx] • (A.9) 

The contribution of single charge exchange (SX) /gx = O'sx/fincl is infered from the mea- 
surements of [103] by simple linear regression, /sx = 0.11 + 6 x 10~^(£'7r/MeV). We fur- 
ther assume a constant ratio of double charge exchange (DX) to scatterings preserving the 
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incident the pion charge (NX), gDx = ctdx/^nx = 0.1. For the multiphcities ^.^^ of the in- 
dividual processes a = SX, NX, and DX we resort to data of inelastic ^He-p scattering. To 
this end, one observes that for incident beam kinetic energies in excess of ~ (50 80) MeV 
(above which also pion inelastic scattering starts dominating over absorpion) the branch- 
ings into ^Henp, T2p, Dn2p and DDp become approximately energy independent with 
respective values ~ 0.4, 0.2, 0.2 and 0.1 (for moderate energies < 300MeV). Adopting 
those values yield the multiplicities for '^He-vr^ reactions as listed in Table 1. Note that 
the ^-values imply a quantitative difference between 7r~-^He and tt^-^Hc scattering. 

In a similar manner as in Section 2.3 where we introduced the correction factor k one 
can account for the incomplete stopping of pions in reactions on ^He by defining an effective 
cross section for tt-^^Hc i , i = ^He, T, D, n and p, averaged over the pion-trajectory, 

= r^M dtPsurv(t)i^7r^4He ^ {'^^±^He,aV) X (A-IO) 
^ Q;=inel,abs 

where 

Psurv(i) = exp (- f ^^Ih^-^^^dA (A.ll) 

is the timc-dilatatcd survival probability of the pion. Note that it was not possible to 
define a correction factor k as in (2.29) because the effective rates for vr^'^He — )■ T and 
TT+^He ^He vanish below the threshold of inelastic 7r='=-^He scattering. 

The output of secondary ^Li, triggered by 7r='=-^He reactions and prior to potential 
destruction on protons, can be tracked via 



dneu 



dt 



"^He E {(C^He^.)*(iVeLi)/. + (r:^lHe^,)*iVeLi} (A.12) 



i=3He,T 



where ^^^^4^^'^^^ = n^j ^'''^l^He^j^^ are the effective rates for energetic i = ^He, T 

production. The dependence of E^^ on time is suggested by the * in (A.12) indicating that 
the effective rates are to be "convolved" with the number of produced secondary ^Li per 
injected A = 3 nucleus, A^bl;. For inelastic scattering this involves an additional average 
over the spectrum /j of primary mass-3 injection energies Ei^\^, 



{Neu)f, = / di^Mn f^{E^,in, E^{t)) / dE, ^ ^^^^/^ ■ (A.13) 

where we employ [14] 

m,^, E.) = y: ^^^i^^'"-'^' X f . (A.14) 



fc=l,2 



Note that we have introduced a kinematical cutoff in the last factor; £'i,max(-E'7r) is the 
maximal kinetic recoil energy of A = 3 in the rest frame of the mother ^He nucleus (which 
approximately coincides with the frame of the thermal bath) and can be found by simple 
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kinematical considerations (see also below). The spectra (A. 14) are obtained from fits to 
experimental p-'^He data with coefficients (i^shc) = (1)0.19), (TsHg/MeV) = (1.9,4.9) and 
(wt) = (1,0.15), (Tt/McV) = (1.3,3.9) [104, 105]. The normaHzation factor J\f is chosen 

such that J^''^'^ dEi;,n fi = l. 

Finally, we remark that for pion absorption, the resulting non-thermal ^Li-yield is 

easier to track because T and ^He are injected mono-energetically in the CM frame so that 
{Nei^i) J, —7- Nei^i as in (2.36). The primary A = 3 injection energy in the rest frame of ^He 
(the frame of the thermal bath) then depends on the scattering angle and lies within the 
interval 



-E'3,m(-£^7r) G 



2m3 



(A.15) 



where /x,r4 and are the respective reduced masses before and after scattering. In our code, 
we adopt the central value. Figure 11 shows the efficiency for production of a secondary 
^Li nucleus, per injected ^He+T as a function of temperature. Pion absorption (solid 
and dotted lines) a 'priori is the most efficient way of generating ^Li as the ^ = 3 recoil 
is maximized. The probabilities have to be weighted by the mass-3 multiplicities and 
relative strengths of absorption and inelastic scattering. Moreover, it has to be stressed that 
Fig. 11 does not include ^Li-destruction on thermal protons which ultimately determines 
the surviving lithium fraction. 

A. 4 Electromagnetic energy injection 

Our treatment of the non-thermal electromagnetic component introduced in the decay of 
pious, kaons, and muons is mainly along the lines of Ref. [106]. The injection of energetic 
primaries and/or 7-particles leads to a quick formation of an electromagnetic cascade 
via the processes of inverse Compton scattering and e^-pair production on background 
photons. The resulting "zeroth-generation" differential photon spectrum is given by a 
broken power law [107], 



Vo(^^/^iow)"^-^ for < Eio 



p^{E^ 



i^o(^7Mow)-'-° iov E,„^ < E^ < Ec (A.16) 
for E^ > Ec 



with a cut-off at the pair production threshold Ec — mg/22T and a power break at 
^low — mg/(80T) [42]. The overall normalization of the spectrum is determined from the 
primary injected energy Eq, 

Eo = J dEyEyPyiE^) = Ko^iow(2 + HEc/Eio^). (A.17) 

The photons in the cascade (A.16) undergo further degradation via the (slower) pro- 
cesses of Compton scattering, pair-production on nuclei, and elastic 7-7 scattering so that 
the total number of energetic photons is given by a competition of injection rate Finj and 
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Figure 11: Number of secondary ^Li produced per injected ^He + T prior to potential destruction 

on protons. The dotted line shows A^en for absorption on threshold [compare (2.36)]. The solid 
and dashed lines show the probability of producing ^Li from in-flight reactions of a pion with 
fixed kinetic energy -B^lint. = 180 MeV on ''He via absorption and inelastic scattering, respectively. 
For the latter compare with (A. 13). These curves have to be convoluted with the actual energy- 
dependent branching into mass-3 elements as in (A. 12). Note that p-burning will have a substantial 
effect on the ultimately surviving fraction of secondary ^Li down to the lowest temperatures shown. 



total energy loss rate r^(£'^). The energy spectrum can then be obtained in form of a 
quasi-static equilibrium solution 



J-y 



nx- 



(A.18) 



where nx is the time-dependent number density of the decaying particles, and Finj = t-^^. 
For the thermalization rate we put 



= Tbh.p + rBH,4He + HE^Wc + 0.5F^^, 



(A.19) 



where Fbh,n, ^c, and F-y-y are the respective rates for Bethe-Heitler scattering, Compton 
scattering, and 7-7 scattering. For Fbh.at = 'riNO'BB.,N we employ the cross section formulae 
as given in [42]. The rate for Compton scattering, Tq = ngCKN, depends on the Klein- 
Nishina cross section ckn [108] and on the total electron (positron) number density Ue, 
with an average fractional energy-loss per scattering [107], 

4/3 



^(^7)-! in(2E^/me) + l/2- 
The rate of photon scattering on blackbody photons, obtained in [109], reads 



2n397r*a^ m, 

3654 2^ \m, 



E-y 

rUe 



(A.20) 



(A.21) 
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and is valid below the e pair production threshold. The factor of 1/2 in (A. 19) can be 
understood by noting that photon-photon scattering tends to split primary gamma rays 
into two gamma rays carrying on average each 50% of the primary energy [109]. Taken 
all together, then defines an approximate survival time for the "zeroth-generation" 
photons. 

The production and destruction of elements due to electromagentic energy injection 
can then be described by the Boltzmann equations for nuclei T,A,P, {At > Aa > Ap): 

poo 

-YaY] / dE^f:^'%E^)a^+A^p{E^). (A.22) 
p Jo 

In our Boltzmann code we include all processes determined in [106] from which we also 
take the photodissociation cross sections a^^i^j . Finally, wc remark that we also take into 

account the secondary production of ^Li by fusion-reactions of energetic mass-3 elements 
which are produced by photo-spallations of ^He. The efficiency of producing ^Li can be 
found in analogy to (2.36). 
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